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EDITORIAL 


PROCEEDINGS OF THE PHYSICAL SOCIETY 
SECTIONS A and B 


In the light of experience gained over the course of the past year it is now possible 
to make a clearer statement of the basis on which division into the two Sections 
of the Proceedings is being made. 

In general, it has been found that the division which best meets the needs of 
Fellows is one which brings papers on microphysics and the physics of elementary 
particles into the one Section, A, and papers on macroscopic physics into the 
other Section, B. While it is essential to maintain a certain amount of flexibility 
in the allocation of papers, a paper will, as far as possible, be put in the Section 
where the main interest lies, e.g. a paper on counter-technique would appear in 
B if the circuitry is the main interest, but would be in A if the applications to actual 


counting are of major importance. 


A list indicating this subject division is given below. 


A, C, STICKLAND, 
Secretary—Editor. 


Section A 


Thermodynamics 
Electrodynamics 
Statistical mechanics 
Quantum theory 
Nuclear physics 
Electron theory : 
dielectrics 
Molecular and atomic structure 
Spectra of atoms and molecules over the whole 
frequency range 


metals, semiconductors, 


Physics of crystals (including luminescence 
as a means of investigating this) 

Photoconductivity 

Theory of solids, liquids and gases 

Magnetism (theoretical) 

Cosmic rays 

Low temperature physics 

Standards relevant to this Section 


Section B 


Thermal, magnetic, electrical, elastic, rheo- 
logical and other mechanical properties of 
matter 

Hydrodynamics and aerodynamics 


Acoustics 

Metallography 

Geometrical optics, optical design and 
microscopy 

Interferometry, diffraction and _ classical 


physical optics generally 

Electric and electronic circuits 

Dielectrics and semiconductors : measurement 
and theory of measurement 

Electric discharges 

Electron optics 

Geophysics 

Radio 

Physics of the atmosphere and ionosphere 

Physics of the sensory processes 
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Colour physics 

Astrophysics and solar physics 

Applications of spectroscopy (e.g. chemical 
analysis, structure analysis) 

Analysis of crystal structure 

Radiography 

Luminescence 

Photoelectricity 

Thermoelectricity 

Piezoelectricity 

Photoconductivity 

Magnetic properties of materials (including 
methods of measurement) 

Production of low temperatures 

High pressure physics 

Vacuum physics 

Thermionic emission 

Standards relevant to this Section 


to 


Heat Transport in Superconductors 


By K. MENDELSSOHN anp J. L. OLSEN 


Clarendon Laboratory, Oxford 


MS. received 6th Fuly 1949 


ABSTRACT. The heat conduction processes in superconductors have been discussed 
on the basis of the analogy with liquid helium II. It is suggested that under certain con- 
ditions the superconductive metal can exhibit a type of heat transport corresponding to the 
heat flow associated with the fountain effect. The heat conductivity of a number of pure 
metals and alloys has been measured in the normal and in the superconductive state and the 
results have been analysed with reference to the hypothesis mentioned above. In addition 
to the change of heat conduction with temperature, the magnetic hysteresis of the heat 
conduction has been investigated. The use of these phenomena as make-and-break 
thermal contacts at very low temperatures has been suggested. 


Ste) UNTRO DUCT ON 
F the various manifestations of superconductivity the phenomena associated 
with the transport of heat appear to be particularly difficult to explain. 
This is due to the fact that the heat resistance of the metal in the super- 
conductive state may be higher ov lower than in the normal state at the same 
temperature. While the former case usually applies to a pure metal and the latter 
to alloys, this in not a general rule. ‘Thus the available experimental material 
presents a somewhat confused aspect which none of the theoretical work has 
attempted to explain. A few years ago one of us (Mendelssohn 1946), when 
discussing the similarity between superconductivity and liquid helium II, 
suggested a possible mechanism which seems capable of accounting for the 
qualitatively different behaviour of pure metals and alloys in heat conduction. 
Indeed, the fact that this difference should have followed as a necessary consequence 
from the analogy between the two phenomena was taken to be strong support 
forts 
In order to test our hypothesis of heat transport in superconductors more 
experimental material than exists at present is required, and it was therefore 
decided to begin systematic experiments on the heat conduction of superconductive 
pure metals and particularly of alloys. A considerable amount of work in this 
field has been carried out in Leiden (cf. de Haas and Bremmer 1936, de Haas and 
Rademakers 1940) in which the dependence on temperature of the heat resistance 
of a number of alloys and pure metals was determined. This was supplemented 
by experiments in this laboratory (Mendelssohn and Pontius 1937) in which the 
change of heat resistance in the transition between the superconductive and the 
normal state was studied in some detail for the alloy Pb 90°,—Bil0°% and for 
pure lead. Finally it was found (Daunt and Mendelssohn 1946a) that the 
temperature gradient in a superconductor is uninfluenced by the setting up of a 
persistent current. This fact, which shows that the Thomson coefficient is zero, 
leads to the important conclusion that the entropy of the electrons in a super- 
current, too, is zero. A survey of this material shows that it would be desirable 
to extend the work on heat resistance to a greater number of substances, and 
especially to combine determinations of the temperature dependence with the 
taking of magnetic cycles similar to the observations of Mendelssohn and Pontius. 
A considerable number of measurements of this kind has now been carried 
out and the present paper gives a survey cf the results obtained so far. In view 
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of the extent of the field of research the work to date is far from complete and can 
only be considered as a first report. However, in view of the conclusions which 
can be drawn already, and because continuation of the work for several years is 
anticipated, it seems desirable to publish the first results at this stage. 


§2. MECHANISMS OF HEAT TRANSPORT 


Our hypothesis of the heat transport phenomena is based on the assumption 
that the electrons in a superconductor can be treated as being made up of two 
inter-penetrating fluids of different entropy (Mendelssohn 1945) similar to 
helium I]. The normal (n-) constituent has ordinary properties, including 
Anite entropy, while the other (z-) constituent has anomalous properties and zero 
entropy. ‘The z-fluid makes its first appearance at the transition point and its 
concentration increases at the expense of the n-fluid as the temperature is lowered. 
Evidence from magnetic determinations (Daunt, Horseman and Mendelssohn 
1939) indicates that at very low temperatures the concentration of the n-fluid in 
a superconductor approaches the value zero. 

Enumerating the processes contributing to the heat transport, we have first 
to consider in the superconductor the effect of the crystal lattice. Unless we care 
to liken it to the conduction of the container, this term has no counterpart in the 
case of helium II. ‘Then comes the conduction of the n-fluid, which is evidently 
quite small in liquid helium but large in the case of electrons. ‘This, of course, 
is the ordinary metallic heat conductivity which ts related to the electrical conduc- 
tion by the Wiedemann—Franz law. Compared with the normal state, this term 
must decrease in the superconductive state with falling temperature as the 
concentration of the n-fluid decreases. It should tend to zero at absolute zero. 
The heat conduction of the z-fluid must be zero as it has zero entropy, but this 
does not mean that the z-fluid cannot contribute to the heat transport. The 
setting up of a temperature difference within the substance is tantamount to 
creating a difference in the concentration of the two fluids. This in turn results 
in a flow of z-fluid towards the hot end*, and in a return flow of n-fluid in the 
opposite direction. A circulation of this type is accompanied by a large heat 
transport. It differs from an ordinary convection by the fact that at the warm 
end the z-fluid is transformed into the n-state and in doing so must take up the 
total entropy at this temperature. This circulation is well known in the case of 
helium II, where it gives rise to the fountain effect (Allen and Jones 1938, 
H. London 1939). It is this process which is responsible for the very high heat 
transport of the liquid. 

The question arises whether a similar circulation can take place in a super- 
conductor. In order to maintain it, a super-current must flow in the direction 
cowards the warm end, and it must be compensated by a normal current in the 
opposite direction. However, here the main difference between helium and a 
superconductor, the electrodynamic properties of the latter, comes into play. 
The super-current must be a surface flow while the induction inside the conductor 
has to be zero. The normal return current and its magnetic field would have to 
be distributed over the whole cross section of the conductor. Because of the 
‘Meissner effect no such field is permitted and the return flow cannot, therefore, 
take place. It is important to realize that a circulation flow similar to that in 

“* The mechanism of this flow of z-fluid has been interpreted by us as a kind of diffusion process 


under its zero point momentum (Mendelssohn 1945), and this assumption has gained support from 


experiments on the helium film (Daunt and Mendelssohn 1946 b). 
I-—2 
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helium with its high heat transport is only prohibited by electrodynamics, and that 
if conditions can be created which allow the induction inside the superconductor 
to be different from zero, this heat transport should take place. 

Such superconductors with finite magnetic induction do in fact exist. The 
systematic investigations carried out in this laboratory on the magnetic and thermal 
properties of superconductors have shown that superconductive alloys are 
electrodynamically inhomogeneous. Their properties can be described as those 
of a magnetic ‘sponge’ (Mendelssohn 1935) whose meshes are made up of super- 
conductive and normal regions. It is, in fact, almost impossible to experiment 
with a superconductor which does not show traces of this sponge structure since 
the magnetic properties appear to be extremely sensitive to chemical and physical 
impurities. According to our views we should therefore expect in impure 
superconductors and in alloys an entirely new heat transport mechanism which 
is exactly analogous to the high heat flow in liquid helium II. The observed. 
fact that the heat transport in pure superconductive lead is smaller than in the 
normal metal, whereas this relation is reversed in lead containing 10°% of bismuth 
(Mendelssohn and Pontius 1937), appears most convincing. 

Thus the heat ‘conduction’ of the metal in the normal and in the super- 
conductive state can be added up as follows. Denoting by L the heat conductivity 
of the lattice, E the ordinary electronic heat conduction, p, the fraction of electrons 
in the superconductive state, and G the heat transport by circulation, similar to 
that in liquid helium IT, we can write for the total heat ‘conductivity’ K of the 
metal : 

(a) in the normal state 


Ku= hak ee eee (1), 
and 
(b) in the superconductive state, for («), B=0 
K=L+i74 =9,), arts 4 
for (B) BAO kK 21 esc ee 3) 
and since /(1—p,) will vary between 0 and 1, 
Ae a eee (4) 
but IG, ee Keg. se) Jeitalsaia) (5) 


In view of the fact that all the terms depend in a complex manner on temperature 
and structure, numerical evaluation of equations (1), (2) and (3) is clearly out of 
the question. However, as equations (4) and (5) show, the expressions for K,, 
and K,seem to lead to a qualitative interpretation of the heat conduction phenomena 
in superconductors. For instance, in a very pure metal where the induction 
is Zero we should always have a smaller heat conduction in the superconductive 
than in the normal state. Moreover, near absolute zero in such a metal the heat 
conduction in the superconductive state must become very small since p, will 
approach the value 1 and the second term in equation (2) disappears. The 
assumption that at absolute zero all conduction electrons in the metal are in the 
superconductive state is based on magnetic determinations of the entropy 
(Daunt, Horseman and Mendelssohn 1939) and appears to be borne out by the 
observed high heat resistance of pure superconductive lead at the lowest temper- 


ature (de Haas and Rademakers 1940). Under these circumstances it is only the 


lattice which contributes to the heat conductivity. 
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Whether in cases where the induction is different from zero K,, or K, is the 
larger must depend, as equation (3) shows, on the relative importance of the 
terms f(1—p,) and G. While a proper assessment cf f(1 —p,) can only be made on 
the basis of a theoretical model of superconductivity, one has reason to believe 
that it should decrease as a monotonic function of the temperature. The ano- 
malous heat transport G will depend on a number of factors. It must be prc- 
portional to Sp,, where S is the total electronic entropy of the metal. Since 
p,—> 9 at the transition point and S-+0 at absolute zero, G must pass through a 
maximum with varying temperature. G is inversely proportional to the 
resistance encountered by the return flow of the circulation in the normal phase, 
but will also depend on the thermo-electrical properties of the latter. 

The whole circulation process might be described as an internal thermocouple 
consisting of a superconductive and a normal lead. No current will be possible, 
for instance, if the thermo-r.M.F. between the warm and the cold end is zero since 
a current in the normal phase depends on the existence of an accelerating 
potential. In other words, under these circumstances G will become zero, and 
we think it possible that this influence of the thermo-electric properties may 
account for the qualitatively different behaviour of different superconductive 
alloys. This conception appears to be borne out by the results on lead—bismuth 
and the lead-tin alloy described in this paper. 


Gia) TRUS SOAS 

The method employed followed the usual pattern of measuring the temperature 
at two places on a rod of the metal along which a steady temperature gradient is 
established. The accuracy at 4°K. was about +2°%, but dropped to +5% at the 
very lowest temperatures. ‘The temperature differences were of the order of 
Q-1°—0-2°. A detailed description of the apparatus employed will be given at 
a later date. 

Tantalum. As an example of a pure metal tantalum was chosen because the 
heat conductivity of this metal has not been measured before. ‘lhe specimen 
used was the same rod of very pure metal whose critical fields had been determined 
by Daunt and Mendelssohn (1937). The results, given in Figure 1, show that 
at all temperatures the heat conduction in the superconductive state is smaller 
than that in the normal state, which is in agreement with equation (4). Moreover, 
the heat conductivity of the superconductive metal decreases rapidly as the 
temperature is lowered, being less than half of the normal at 2:75°K. Also, the 
curve for the normal state appears to be the continuous one, from which the super- 
conductive curve breaks away at the transition point. All these features agree 
with previous work on pure metals as well as with our equations (1) and (2). 
Near the transition point the scatter of our values is considerable and may be 
larger than is warranted by the experimental errors. A possible explanation of 
this will be discussed below. 

Pb 90%-Bi 10%. This alloy, the results for which are given in Figure 2, 
was of identical composition with the sample used by Mendelssohn and Pontius, 
As in the pure metal, the heat conductivity in the normal state follows a continuous 
curve. However, the heat conductivity of the superconductive metal is entirely 

different. Instead of falling below the normal curve it rises steeply. ‘The whole 
pattern is so different from the pure metals that it appears clear that a new type of 
heat transport must be operative. At 5°x. the superconductive heat conductivity 
is three times as large as the normal. Because of the high threshold fields of this 
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alloy the normal curve could not be determined to lower temperatures, but extra- 
polation suggests that at 3°K. the superconductive value for the heat flow may be 
six times higher than the normal one. In this temperature region the super- 
conductive curve shows a maximum, indicating a rapid drop to lower temperatures. 
Since the normal curve is not known, the question has to be left open whether it 
will also show a drop. However, it may well be that the observed maximum 1s 
that predicted in the previous section for the circulation flow, being due to the 
drop in the electronic entropy. 
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Figure 1. Heat conductivity of tantalum. The continuous curve denotes the normal state of the 
metal, the broken curve the superconducted state. © in zero magnetic field; @ in fields 
exceeding the threshold value. 
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Figure 2. Heat conductivity of Pb—Bi 10% alloy. 


Pb 70°%—Sn 30%. In order to test our assumption concerning the influence 
of the thermc-electric properties on the circulation flow a lead—tin alloy was 
measured because it is known (Keesom and Matthijs 1938) that the thermo- 
E.M.F. between these metals in the helium region is very small. Although the 
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therm:-E.M.F, in the alloy is likely to be different from that between the pure 
constituents, one may expect that it is much smaller than between lead and 
bismuth. The results, shown in Figure 3, indeed seem to support this assumption. 
The superconductive curve lies at all temperatures below the normal, and the 
results altogether resemble much more those for tantalum than those obtained on 
the lead—bismuth alloy. In particular, we may assume that the induction in this 
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Figure 3. Heat conductivity of Pb-Sn 30% alloy. 


case was different from zero, and there should, therefore, be a chance for the 
establishment of anormal current inthe metal. That inspite of this no appreciable 
heat flow by circulation takes place can be taken as a strong indication for the 
expected thermo-electric influence. 

Columbium. It seemed interesting to investigate this metal not only because 
of its similarity to tantalum but also because it has a very high transition point 
(9-2°x.). The rod used was of high purity (‘H.S.’ Lab. No. 1919). The results 
(Figure 4) are complex and differ from all previous observations in that the curves 
cross over at low temperatures. Immediately below the transition point the 
superconductive metal, similar to tantalum, exhibits a lower heat conduction than 
the normal. However, as the temperature is lowered the difference between the 
two curves decreases again, and below 3-8°k. the superconductive metal becomes 
the better conductor of heat. It would seem very difficult to account for this 
peculiar behaviour except by two different, competing processes in the heat flow 
mechanism. We have seen that, according to our hypothesis, such a mechanism 
actually exists in the relative influence of the terms f(1 —p,) and G in equation (3). 
Inavery pure metal the latter term should always be zero, but it is well known that 
the induction in the ‘hard’ superconductors is not always zero. In fact it has 
been found that a similar columbium rod exhibited more and more the magnetic 
properties of an alloy as the temperature was lowered (Daunt and Mendelssohn 
1937). We had therefore reason to suspect that the crossing over of the curves 
might be due to the appearance of circulation heat flow at low temperatures. 
In order to investigate this question more thoroughly, the dependence of the heat 
flow on the magnetic field was measured in some detail. 


8 K. Mendelssohn and Ff. L. Olsen 


The magnetic cycles, given in Figure 5, represent the change of thermal 
resistance (in arbitrary units) with transverse magneiic field. The dotted curves 
which give the hysteresis cycles after superconductivity had been destroyed are 
similar at the two different temperatures (4-4° kK. and 2-7°K.), except that there is 
somewhat less hysteresis at the higher temperatures. This is exactly what would 
be expected on the strength of the magnetic measurements. It shows that the 
metal behaves to some extent like an alloy and that this effect increases as the 
temperature is lowered. The striking feature cf these observations, however, 
is the full curve denoting the first destruction of superconductivity. At 44°k., 
where the superconductive heat flow is still lower than the normal, we might 
expect a hysteresis figure similar to that observed by Mendelssohn and Pontius 
on Pb-Bi except that it should be completely inverted. That means in our 
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Figure 4. Heat conductivity of columbium. 


present figure the beginning of the virgin curve should lie somewhat above the 
two peaks of the hysteresis loop. Instead, however, the virgin curve starts well 
below this value, which means that during the hysteresis the heat resistance 
(at about 1,000 gauss) is higher than in either the superconductive or the normal 
state. ‘T’his feature is even more pronounced in the hysteresis figure at 2:7°K., 
where the beginning of the virgin curve is depressed below any other value in 
the compiete cycle.* ‘These results seem to show without doubt that the crossing 
over of the curves in Figure 4 is due to changes of.the heat flow mechanism in 
the superconductive and not in the normal state. The two hysteresis figures 
demonstrate the appearance of some additional form of heat transport as the 
magnetic behaviour of the metal alters with temperature. The evidence offered 
thus appears to support the hypothetical circulation process. 

Pb 99-9°%,-Bi0-1%. Asa further test of our hypothesis it appeared tempting 
to try to create artificially similar conditions to those which according to our 
analysis obtained in the columbium specimen. In order to simulate these 


* It is possible from this figure to extrapolate the heat conductivity in the superconductive state 


without the circulation process, i.e. to place the start of the virgin curve above the maxima. This 
value is marked by a cross in Figure 4. 
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Figure 5. Hysteresis figure (heat resistance W in arbitrary units against magnetic field) 
for columbium. 
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Figure 6. Heat conductivity of Pb-Bi 0°1% alloy. 
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IO 
conditions we decided to add a very small amount of a second constituent to a pure 
metal. For this purpose we chose lead, since its heat conductivity had been 
measured in Leiden and in Oxford, and the results indicated a close approximation 


to ideal conditions. To this 0-1°% cf bismuth was added. The results (Figure 6) 


show indeed that, similar to our columbium rod, the difference between the 


superconductive and the normal curve first increases with lowering temperature 
In this case the curves do not actually intersect, the: 


and then declines again. 
superconductive one being lower all the time. However, reasonable extra- 
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Figure 7. Hysteresis figure for Pb—Bi 0:1% alloy. 


polation suggests that the curves may possibly cross over at about 2° k. 


curves at a higher temperature. 


A slightly 


higher percentage of bismuth would most likely produce an intersection of the- 
An even more striking similarity with the 


observations on columbium is provided by the hysteresis curves (Figure 7). 
When comparing the curves denoting the first destruction of superconductivity 
it can be seen that again, as the temperature is lowered, the beginning of this 
curve is depressed with respect to the rest of the hysteresis cycle. It appears. 
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almost certain that this anomalous rise in the superconductive heat conduction 
must be the first indication of the circulation process which in the alloy with 10%, 
bismuth completely overshadows all other effects. 


S44; CONCLUS ON 

As already stated, the work has to be regarded as being of a preliminary character, 
and the results obtained so far emphasize the need for further and more detailed 
information. On the other hand, the available data seem to provide fairly 
satisfactory evidence for our qualitative model of heat transport phenomena,. 
especially the existence of a new type of heat transport which can under certain 
(controllable) circumstances take place in the superconductive metal appears 
well established. Whether this process is in fact identical with the circulation 
flow postulated by us and to what extent it is the counterpart to the high heat 
conduction in liquid helium II can only be decided with certainty when more 
data are available. It should be pointed out, however, that all the results obtained 
so far can be regarded as confirmatory. In particular, the fact that the complex 
and unexpected results on columbium have proved amenable to our interpreta- 
tion is most encouraging. ‘The positive outcome of the two test experiments, 
one (Pb—Bi 0-1°,) simulating the conditions in columbium and the other (Pb-Sn) 
the influence of the thermo-electric effects, must also be counted as supporting 
evidence. However, further experiments on these lines are clearly required. 

In order to summarize our results we have plotted the ratio of the super-— 
conductive to the normal heat conduction K,/K,, against the reduced temperature, 
T/T., where 7, is the transition temperature. In this way our results can be 
compared with the theory of Heisenberg (1948). The deviation of our curves. 
from the theoretical one can be seen in Figure 8. Apart from tantalum, all our 
results indicate higher ratios, K,/K,, than the theory predicts. In view of our 
hypothesis, this is not surprising. Heisenberg’s theory only considers the 
process corresponding to our equation (2), making no provision for the circulation 
transport postulated by us. It is therefore interesting to note that all our ratios 
are higher except for tantalum where, according to our conceptions, the term 
Gshould be zero. Moreover, Pb—Sn, which should have the smallest G, also shows 
the smallest excess, while for Pb—Bi 10°%, which will have the highest G, the excess 
ratio is so large that the scale had to be changed in order to accommodate its curve 
in the same diagram. These curves, too, therefore seem to provide evidence for 
a new process carrying large amounts of heat in superconductors. Quite recently 
Hulm (1949) has summarized some heat conductivity measurements in this type 
of diagram and concluded that the ratio K,/K, may be a characteristic function 
of temperature. In view of our results and of those obtained by Mendelssohn 
and Pontius we cannot but regard the mutual agreement between Hulm’s curves 
as rather accidental. It would seem that in none of the cases investigated by him 
was G appreciable, but an analysis cf his results suggests that, if he had further 
increased the mercury percentage of his Sn—Hg alloys, he might possibly have 
obtained a ratio K,/K, larger than unity. 

Finally two possible future developments of this work should be mentioned. 
Firstly, in our treatment of the term G we have omitted the possibility that the 
circulation process may be limited at a given temperature by the superconductive 
threshold current. This would result in a saturation effect and the heat 
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‘conduction’ should be larger for small heat currents. No such effect has as 
yet been established, but it may be responsible for the scatter in the tantalum results 
mentioned above as well as for the unexplained scatter observed by de Haas and 
Rademakers (1940) in lead. Secondly, the phenomena of heat flow in super- 
‘conductors may find a practical application in experiments at very low temperatures. 
Any cooling process carried out in successive stages will require make-and-break 
thermal contacts. On the basis of their recent experiments Hudson, Hunt and 
Kurti (1949) have suggested for this purpose the use of a thin layer of liquid 
helium between metal contacts. The results described in the present paper now 
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Figure 8. Ratio of heat conductivities against reduced temperature. 


seem to offer an alternative or possibly a supplementary solution to this difficult 
problem. The difference between the heat conduction of a metal in the super- 
‘conductive and the normal state can in certain cases be considerable, and one has 
reason to expect that, for instance in pure metals, at very low temperatures the 
ratio may be of the order of one to twenty or more. Furthermore, judicious 
-adjustment of the amount of asecond constituent can probably produce temperature 
variation of the ratio in such. a manner asttoneurt the particular experiment. 
-Make-and-break thermal contacts consisting of rigid connections of a metal whose 
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heat conduction can be altered by a magnetic field may therefore often be found 
to be sufficient. ‘They would have the advantage of being free of moving parts, 
and by a suitable choice of the metal the necessary magnetizing field could in many 
cases be kept small. ; 
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ABSTRACT. The magnetization curves of thin superconducting tin plates were measured: 
in transverse magnetic fields in order to investigate the nature of the intermediate state in 
such specimens. The curves showed peak magnetizations much higher than H,/47, indi- 
cating the presence of fields much greater than the critical field at the edge of the plates.. 
This effect is qualitatively similar to that predicted by Landau’s theory of the intermediate 
state, but quantitative agreement is not obtained. Resistance measurements on thin strips 
of tin in transverse fields showed that resistance only reappears for fields considerably 
greater than those required to start the destruction of superconductivity, 


S15 INTRODUCTION 

HE experiments of Désirant and Shoenberg (1948) on the intermediate 

state of thin superconducting tin and mercury cylinders in transverse 

magnetic fields indicate qualitative, and to some extent, quantitative 
agreement with the model of the intermediate state first suggested by Landau 
(1943), and worked out in more detail by Andrew (1948), Certain features of the 
theory suggested, however, that it would be desirable to extend the experinients to. 
superconducting plates in transverse fields in order to test the theory further. 


* Now at Department of Natural Philosophy, University of St. Andrews, Scotland. 
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For an ellipsoid whose dimensions are large enough to allow neglect cf all 
surface energy effects, simple argument shows that the magnetization curve 
consists of two straight lines. For 0<h<H,(1—n), where h is the applied 
magnetic field, H, the critical field and 47” the demagnetizing coefficient, the 
ellipsoid is in the superconducting state, and the intensity of magnetization J is 


given by 


—4rl=h/(V=n). =) eee (1) 
For H,(1—n) <h<H, the ellipsoid is in the intermediate state, and / is given by 
45 )=(Ho hii 9 Se ea eee (2) 


For h>H, the ellipsoid is in the normal state, and /=0. ‘This magnetization 
curve will be described as the ‘bulk curve’. For the very oblate ellipsoids to which 
the flat plates used in the experiments approximate, the value of 7 is nearly unity 
and the initial slope of the ‘bulk curve’ is very great, corresponding to the large 
value of 1/(1—7). 

The magnetization curve for an ellipsoid for which surface effects are not 
negligible, differs from the bulk curve in the manner indicated in Figure 7, where 
the bulk curve is shown by a broken line. The salient features predicted by the 
theory, which will be discussed in more detail in § 3, are: (1) the considerable delay 
in formation of the intermediate state, giving rise to a pronounced peak; (1i) a 
steeper slope for the intermediate state portion of the curve than that for the bulk 
curve; (iii) the attainment of the normal state when the applied field 1s appreciably 
less than the critical field. 

The magnetization curves of thin tin plates of various thicknesses have been 
_ determined at various temperatures, and since it was possible to use specimens in 
which the range of thicknesses extended to much smaller values than was the case 
with the cylinders used in previous experiments, the departure of the curves from 
the ‘bulk’ shape was expected to be considerably greater. This was in fact found 
and, as with cylinders, qualitative agreement with the theory in certain respects was 
obtained ; in other respects, however, more serious disagreement was found. 


§2; EXPERIMENTAL DETAILS 


All the plates were prepared by rolling spectroscopically pure tin (Johnson and 
Matthey 2356) to the required thickness. Discs were then cut out, usually with a 
clean cork borer. One specimen, Sp. 12, was of oblong shape, and this was a 
small portion of a specimen E3 used in the resistance measurements described 
later, for which the ratio of the resistance just above the transition temperature 
(3-8°k.) to that at 18°c. was found to be 3-0 x 10-4, indicating a high degree of 
purity. The discs were about 8 mm. in diameter and their thickness varied from 
2 to 400 p. 

In order to prepare the thinnest specimens it was found convenient to allow 
the strip of tin which was being rolled to stick lightly to the rollers as it passed 
through, as in this way the smoothest surface could be obtained; most of the tiny 
wrinkles which appeared could be removed, after the foil had been rolled to the 
required thickness, by widening the gap between the rollers very slightly and 
passing the foil through once more. 

In all cases the mean thickness of the specimen was calculated from measure- 
ments of weight and surface area. 
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The method of measuring the magnetic moments of the specimens differed in 
only one way from that used by Désirant and Shoenberg (1948). In the present 
experiments the initial slope of the magnetization curve was generally so great 
that the current required for the compensating coil was greater than the solenoid 
current itself, so that it could not, as in the previous experiments, be taken as a 
fraction of the solenoid current, but was supplied instead from a separate circuit 
and was measured directly on a milliammeter. The overall sensitivity of the 
arrangement was of the order of 3 x 10° mm. scale deflection per £.M.U. of magnetic 
moment. For the thickest specimens the photoelectric amplifier and second 
galvanometer were dispensed with, and the magnetic moment was measured 
directly as a scale deflection of the first galvanometer, using the compensating coil 
merely for calibration. 

All temperatures quoted are deduced from the vapour pressure of the helium 
bath by means of the 1932 Leiden scale. 

The mounting of the foils used for the resistance measurements is described 
elsewhere (Andrew 1949). 


§3. MAGNETIZATION CURVES 
3.1. Experimental Results 


Magnetization curves were measured at temperatures of 3-49°k. and 2:99°K, 
for all specimens, and at 2:50°K. for some specimens. ‘Typical results are shown 
in Figures 1—6, in which the ordinates are 47 times the magnetization in gauss and 
the abscissae the applied field in gauss; the full curves are calculated from the 
theory given below, using the values of the parameter A’ indicated on p. 22. It 
will be seen that all the magnetization curves show a considerable peak for low 
values of the applied field. Certain quantitative features of the curves are 
exhibited in Table 1 for comparison with theory. 

The critical fields of specimens Sp. 5, 12, 15, and 16 have been measured at a 
number of temperatures in magnetic fields parallel to their surfaces, and from these 
measurements it is found that the deviations of the transition temperatures of these 
specimens from 3-72°K. are not greater than 0-01°k. in any case. The critical 
fields of the other specimens have not been measured, since they are not of decisive 
importance in comparison of the results with the theory. 


3.2. Theory 


A theory of the intermediate state has been worked out by Landau (1943) and 
Andrew (1948), and therefore only the results needed here will be quoted. Bya 
calculation similar to that for transverse cylinders, it is found that the point of 
entry into the intermediate state from the superconducting state should take place 
at an applied field pH, where p is greater than the value (1 —) for the bulk curve. 
‘The values of p for the laminar and thread models are given by the following 


relations: 
Laminar model py=v+2(4Nv)*?, = «22... (3) 


Whreadimodel = ppv —2(Mey, wanes (4) 


where v is written for (1—#), and 
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Figure1. Sp. 11,2-988°x., A’/L=1-72x10-1.  Figure2. Sp. 16, 2:984°x., A/L=1-:0x 10-3. 
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Figure 3. Sp. 15, 2:996° x., A’//L=4-2 x 107%. 
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Figure 4. Sp. 15, 3-492°x., A’/L=7-1 x 10-3. 
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Figure 5. Sp. 5, 3:490° x., A’/L=1:55 x 10-1, Figure 6. Sp. 12, 3-514°x., A’/L=8:-0x 10-3. 


Figures 1-6. ‘Transverse magnetization curves. In each case the full curve is theoretical, assuming 


values of A’/L indicated. ~-—+- discontinuous rise of —4aI; © increasing fields ; 
+ decreasing fields. 
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A‘ \2'5 

Thus py! =1+3-1 (=) hee ae ER eee (5) 
A’ \12 

and priv=l 45-9 (5) ‘Stee ee Oa vector eters (6) 


The complete magnetic behaviour of an ellipsoid can be derived from the 
three equations 


Hoahetinlt “Beatle 60 os ee (7) 


For an ellipsoid in the intermediate state, ®, the free energy multiplied by 47 may 
be expressed as H,.B +¢, where ¢ is a term which takes interphase surface effects 
into account. For the laminar model described by Landau, 


_ 3BYH.(H.— B)?? A’ 2/3 
242 Ib, ; 


where L is the thickness of the material in the direction of the applied field, and 
A’ is given by 
S70 


H 2 


c 


N= 


= Nn Sh Be ae ees (9) 
where A is the penetration depth of the field into a superconducting domain, and 


x is the surface energy per unit area between the normal and superconducting 
phases. A similar expression for ¢ is derived from the thread model described by 


Andrew. 
“Ht 
-401 


Belt 


h/F, 


Figure 7. Theoretical magnetization curve. Full curve for INE = NO, p= (02 
Broken curve for A’/L=0 (‘ bulk curve’). 


In applying the theory to long cylindrical specimens perpendicular to the field, 
721s put equal to 3, but it is difficult to decide what value to take for L in the formula 
for ¢ given above, since it varies from 27 in the centre of the specimen to zero at its 
edges, where 7 is the radius of the cylinder. The theory can, however, be applied 
directly to the similar problem of thin plates perpendicular to the field, as here no 
doubt arises as to the value of L to beinserted. Another difficulty arises, however, 
because a uniform flat plate with sharp edges does not correspond exactly to the 
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pertect elltpsoid required to give uniform magnetization, so that its behaviour 
cannot be described completely in terms of a simple demagnetizing coefficient and, 
as will be mentioned later, the experimental values of the initial gradients cf the 
magnetization curves do in fact differ from those calculated on the assumption that 
the uniform plate may be treated as its inscribed ellipsoid. 

In the limit A’/L 0, equations (7) and (8) lead to the bulk curve already 
mentioned in §1. For A’/L =10-* and v=10-, these equations lead to the curve 
shown in Figure 7 ; two discontinuous jumps are predicted at fields hk, andh,. The 
first of these occurs when the free energy of the superconducting state becomes 
just equal to that of the intermediate state, so that h,=pH,, where p is given by 
equation (8). The second jump AB corresponds to the fact that at a field h, the 
normal state becomes energetically favourable but its attainment would require 
the sudden expulsion of the remaining superconducting domains in the specimen. 
Experimentally this is found not to occur, and the curves agree much better with 
the assumption that the intermediate state persists in a metastable condition giving 
rise to the continuation of the magnetization curve along AC. 


3.3. Comparison of Theory and Experiment 


From (1) the initial gradient of the magnetization curve, —47//h, should have 
the value 1/(1—m7). The demagnetizing coefficient of an oblate ellipsoid of 
revolution of diameter d and thickness L is given by 


ap IL, L\2 LN8 
Peo ee ee a 


"This relation may be used to calculate the expected value of the initial gradient if it 
is assumed that the discs behave like thetr inscribed ellipsoids. ‘The experi- 
mentally determined slope was, however, always about 2/3 of the calculated slope 
for the discs. The reason for this discrepancy is probably as follows. Consider 
a superconducting specimen in the form of an oblate ellipsoid of revolution thin 
enough for its demagnetizing coefficient in a field parallel to its axis to be given by 
the relation 


‘where L is the length of its minor axis and d that of the major axes. ‘Then in an 
applied magnetic field h its magnetic moment M will be given by 

1 Sond OLN PTT IS a8 tages. 

—4aM=305(5) ah pa eee (11) 
‘Thus the magnetic moment of such a specimen in a given field is independent of 
its thickness; physically this means that the distortion of the magnetic field due to 
the presence of such a superconducting specimen depends only on its dimensions 
in the plane perpendicular to the applied magnetic field and not on its dimensions 
in the direction of that field. It seems reasonable to suppose that for a thin 
superconducting plate this result holds even if the cross section of the specimen 1S 
not elliptical, and in particular that it holds for a thin cylindrical disc; if this 1s so, 
then for such a disc the magnetic moment ™ in a given field will be the same as 
‘that of its inscribed ellipsoid. 
i a 
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Since the volume of the disc is 17d? the slope of the magnetization curve will 
be 


eee (0) 


—4nM 4d. 2 oA 
ind@ih 37L 31-n’ 
which is just the value found experimentally. In analysing the results it is 
this experimental value of the demagnetizing coefficient which has been used. 
throughout. 

It will be seen that except for the thickest specimens, after an initial steep rise 
at constant gradient the peak rounds off somewhat, so that the summit is not 
attained until a field has been reached considerably higher than would have been 
required had the gradient remained constant. This effect is on the whole more 
pronounced in the thinnest specimens, and is presumably due to field penetration 
at the sharp edges. In comparing the results with the theory, the ratios 
—4rTnax/H, where Iya, is the maximum magnetization, have been compared 
with the theoretical values of p/v, which is of course equal to —47J,,,x/H, for the 
theoretical curve. 


"7-0 -60 -50 -4:0 -3:0 -2-0 


Figure 8. Variation of (—47Imnax/He—1) with Lv plotted logarithmically. 
Full curve is theoretical, assuming A’=5 x 10~° cm. ; 
LX 225K ©) 299 ace OA Oe 


The extent to which the experimental values of —47/,,,./H, agree with the 
predictions for the laminar model is shown by plotting log[—47/,,,,/H,—1] 
against log Lv in Figure 8. According to equation (5), for the laminar theory such 
a plot should yield a straight line with gradient — 2; such a straight line, taking 
A’=5 x 10-°cm., is drawn on the graph for comparison and it is evident that the 
functional relation is markedly different from that predicted, and the disagreement 
with theory becomes particularly pronounced for very large values of L, that is 
to say for thick discs, if we keep the diameter constant. It is this discrepancy 
which is disturbing from the point of view of the theory, since it is just these 
specimens which yield the most reliable and the most reproducible experimental 
results; as can be seen from the magnetization curve of Sp. 16 at 2:984°x., there 
is no appreciable penetration of the field right up to the top of the sharp peak. 
The thread model gives no better agreement. 

To decide whether the exact shape of the edges of the specimen was critical, 


two discs, Sp. 15 and Sp. 16, were polished electrolytically ; in this way specimens 
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with considerably rounded edges and smooth surfaces were obtained, though 
slight electrolytic etching as well as electrolytic polishing took place. The 
magnetization curves obtained, however, indicated no significant change either in 
the initial gradient of the curve or in the value of —47/,,,,/H,, as compared with 
unpolished specimens of the same thickness and diameter. 

Another feature of the experimental magnetization curves which can be 
related quantitatively to the theory concerns the slope of the falling portion of the 
magnetization curve. From the equations (7) the gradient d//dh of the magnetiz- 
ation curve is given by 


dB apy 
4a(1—n) 2 +1= 5, =| 40-9 55 fae Sp geen avs (13) 


For the laminar model 


ds 21/3 A‘\2/3 - 
reesei ala) ee we 
and for the thread model 
3B 
ea | 4 
dh 27 Par (1 7) A’\28 (15) 
dB? 3(4/2—1)28ya8(1— BLE J 


where y= B/H, and 8 =(2B/7H,)?. For B=4H, both these expressions become 
very nearly d?¢/dB? = —6-1(A’/L)?*. 
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Figure 9. To illustrate how o/(1—n) is obtained. 


If we draw the tangent to the magnetization curve at the point P (Figure 9) 
given by B=4H, and produce it to cut the initial straight line portion of the curve 
at a field h =hz and the axis of field at fg, and define o as the ratio of hy/ha, then 


o=(1—n)| 1- oe | (i E +61 (z)"]. ae (16) 


Although o and (1 —n) are both small, their ratio o/(1 —n) can be calculated by 
a simple geometrical argument (see Figure 9). We have o=ON/OQ and 
1—n=ON/RN;; hence 
itd = IRNIOOK. ag. 8 o/s epider (17) 
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The relative scale along the axes OX and OY for an experimental magnetization. 
curve is independent of the exact critical field of the specimen and depends only 
on the mass of the specimen, the solenoid constant, and the absolute calibration of 
magnetic moment, all of which are known to within trae the values of 
o/(1—n)—1 calculated from the experimental curves have been plotted against 
L-23 for temperatures of 3-49 and 2-99°x. in Figure 10. The points lie fairly 
accurately on straight lines passing through the origin, and from the gradients of 
these lines the values of A’ obtained are 6-3 x 10->cm. and 3-7x10-°cm. at 
3-49° x. and 2-99° k. respectively. These values are in good agreement with those 
obtained by the same procedure as applied to transverse cylinders (Desirant and. 
Shoenberg 1948), namely 6:2 x 10-*cm. and 3-6 x 10->cm. respectively. The 
theoretical magnetization curves obtained using these values of A’ are shown in 
Figures 1-6, and over the intermediate state range where the magnetization has. 
fallen from its high value at the initial peak the agreement is seen to be fairly good. 
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Figure 10. Variation of o/(1—n)—1 with L~*. 


We have so far described only that portion of the magnetization curve traced 
out as the applied field is raised monotonically from zero to a value greater than the 
critical field. On lowering the field to zero again considerable hysteresis is found ; 
for most specimens the field can be reduced to a value lower than that at which 
the rising field curve cuts the axis before magnetization reappears, and when it 
does so it takes place in a jump indicated by a sudden kick of the galvanometer. 
This ‘super-cooling’ effect is similar to that observed in the earlier experiments on 
transverse cylinders. On further reduction of the magnetic field the magnetization 
curve retraces quite accurately the rising field curve until the base of the initial 
peak is reached, when it begins to fall fairly sharply and drops to a negative value 
when the field becomes zero. This ‘frozen-in moment’ is largest for the thinnest 
specimens, for some of which it reaches a value whose magnitude is actually 
greater than that obtaining at the peak of the rising field curve. It is presumably 
associated with the formation of a stable superconducting ring at the edge of the 
disc in which a permanent current can flow, and it is noticeable that this hysteresis 
effect is in fact very similar to that observed with a superconducting annulus. The 
polished specimen Sp. 16 had a smaller frozen-in moment than the unpolished 
specimen Sp. 17 of about the same thickness, and this suggests that the formation 
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of a superconducting ring is favoured by the presence of sharp corners, in agree- 
ment with the results obtained by Shoenberg (1937) for cylinders whose length 
was about the same as their diameter. It is also in agreement with the much 
higher frozen-in moments obtained with the thinnest disc specimens. 

If the formation of the intermediate state took place reversibly over the whole 
range of applied fields, it can be shown quite generally that the area under the 
(J, h) curve should be exactly H,?/87 regardless of the particular assumptions which 
are made concerning the structure of the intermediate state. The theoretical 
curves have in fact just this area, but it is immediately evident that the experimental 
results yield a magnetization curve with a far greater area, the discrepancy increas- 
ing as we go to thinner specimens. Thus for Sp. 11 at 2:98°x. the area of the 
increasing field curve is 7:0 x H.?/87. Now irreversible changes will give rise to 
just such an increase in the area under the magnetization curve as is observed, and 
since most of the difference is associated with the large initial peak of the curve, and 


RR normal 
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Figure 11. Resistance-field curves for E3 at 3:510° kK., using various measuring currents. 
eset icin, 


since the drop from this peak is much more gradual than the theory predicts, we 
must assume that over this range of fields the transitions taking place are not 
reversible. No theory based on a structure of the intermediate state in which the 
changes taking place are reversible will therefore be adequate to describe the 
observed behaviour in this part of the magnetization curve, which must instead be 
ascribed to the formation of some kind of metastable structure of normal and 
superconducting domains which breaks up into a stable ‘intermediate state’ at 
higher fields. The fact that the decreasing field curve follows the rising field 
curve accurately for such fields proves that here we are dealing with reversible 
changes and a stable structure. 


§4. RESISTANCE MEASUREMENTS 
Resistance measurements were made on a number of transverse tin foils in the 
form of narrow strips about 3 cm. long and 1 mm. wide, varying in thickness from 
3-4.to78p. Resistance-field curves for the 78 p foil, of which a portion was used 
as Sp. 12 in a magnetization experiment (Figure 6), are shown in igure 11 for the 
same temperature and are seen to depend strongly on the measuring current. For 
a weak measuring current the foil had a superconducting path for applied fields up 
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to 0:63H,, a value considerably greater than that at which superconductivity begins 
to be destroyed in Sp. 12 as shown by the peak of the magnetization curve; a 
similar high value of the field required to cause the reappearance of resistance was 
found for all the foils. Stronger measuring currents displace the point at which 
resistance first appears to lower fields. : 

Since the magnetization measurements were made in the absence of a measuring 
current they should be compared with the resistance measurements using the 
smallest current, for which the corresponding field at the edges of the foil was only 
about 0-01H,. The absence of resistance in a field greater than that required to 
start the destruction of superconductivity may be due to the normal domains 
growing from the edges and running in from both sides of the strip; in this way a 
central superconducting path can exist along the whole length of the foil until a 
normal domain actually cuts across its entire cross section. ‘Thus no evidence can 
be obtained about the initial destruction of superconductivity in plates from such 
resistance-field measurements. On the other hand, the restoration of complete 
normal resistance seems to occur as expected at a field very nearly equal to that 
required to destroy the last traces of superconductivity as indicated by the magnetiz- 
ation reaching zero. Thus in Figure 6 the point M marks the value of the field 
required to restore normal resistance in the foil of the same thickness from which 
Sp. 12 was made. 


§5. DISCUSSION 


The theory of the intermediate state considered above is based fundamentally 
on the hypothesis that the macroscopic laminae or threads of normal material 
undergo progressive branching until at the surface there results an exceedingly 
finely divided mixed phase. Recent experiments by Meshkovsky and Shalnikov 
(1947 b), in which the field distribution over the surface of a sphere in the inter- 
mediate state was investigated by means of a very fine bismuth probe, show instead 
irregularities in the field on just the same macroscopic scale as were found in their 
previous investigation (Meshkovsky and Shalnikov 1947 a) of the field inside a gap 
cut across a diametrical plane of the sphere. ‘This result appears to disprove the 
branching hypothesis for such a specimen, and it is certainly not more likely to be 
true for the thinnest specimens used in the present experiments, in view of the 
small distance available for such branching to take place. If this is so, the theory of 
Landau and Andrew must at least require considerable modification. Where 
agreement between this theory and the results of the present experiments is found, 
and particularly in the slopes of the magnetization curves in the nearly linear 
region shortly before the magnetization has fallen to zero, it is possible that the 
role played in the theory by the branching of the laminae is not decisive. In this 
case the close agreement between the values of A’ derived from the o values obtained 
in the present experiments and from those obtained by Désirant and Shoenberg 
in the experiments on cylinders may indicate that in the correct theory of the 
intermediate state the surface energy between the normal and superconducting 
phases will play a part similar to that which it plays in the Landau theory. The 
possibility that the agreement mentioned above is fortuitous must not, however, be 
neglected. 

As already pointed out, the large area of the magnetization curves of the thin tin 
plates indicates irreversibility, especially in the region where the magnetization 
begins to fall rapidly. Although the area discrepancy is much smaller in the 
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experiments of Désirant and Shoenberg on transverse cylinders, it is probable 
that the ‘horns’ in the magnetization curves found by them are manifestations on a 
smaller scale of the same process which gives rise to the large peaks in the magnetiz- 
ation curves obtained in the present experiments. Thus little reliance can be 
placed on the values of A’ derived from the p values for thin cylinders in terms of the 
equilibrium theory, and perhaps the higher A’ values suggested by the present 
work are nearer the truth. 
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ABSTRACT. A modified perturbation technique is described for problems in which 
second-order effects are comparable in magnitude with first-order effects, where orthodox 
methods break down. It is applied to the energy levels of paramagnetic ions in a crystal, 
giving an effective Hamiltonian in which the Stark splitting, the anomalous g-value and the 
temperature-independent paramagnetism are clearly exhibited. 


Sele NERO DIL COIN 

N some applications of perturbation theory it happens that first and second 

order effects are comparable in magnitude. The straightforward pertur- 

bation development then leads to difficulties if the unperturbed state is 
degenerate, for the first step is to set up the first-order secular equation, from 
which are derived the zero-order states: while actually a second-order secular 
equation is necessary to get them correctly in this case. 

Such a state of affairs arises in the theory of paramagnetic crystals. Here one 
has to study the energy levels of paramagnetic ions under the influence of the 
crystalline electric field, the spin-orbit coupling and the external magnetic field. 
‘The problem is usefully attacked in stages, first studying the energy levels under 
the influence of electrostatic fields only, and ignoring spin forces; and then, 
with this as the unperturbed system, treating the spin-orbit coupling and external 
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field as perturbations. The contribution of the spin-orbit coupling vanishes. 
in first order, if the unperturbed orbital level is non-degenerate, leaving only the: 
contribution from the external field. This may, in practical cases, according 
to the strength of the field, be larger than, comparable with or smaller than the 
second-order contribution from the spin-orbit coupling. Higher order terms. 
are smaller. 

The purpose of this note is to show how to set up the required second-order 
secular equation in a form applicable to the paramagnetic problem. Here the 
unperturbed system is characterized by certain configurational variables and has 
energy levels which may conveniently be called orbital levels. Owing to the spin 
each orbital level is (2.8 + 1)-fold degenerate. It will be assumed that the lowest 
orbital level, in which one is primarily interested, is non-degenerate apart from the 
spin degeneracy. The perturbation Hamiltonian is 


HLS 8(42S8), eee (1) 


It will be shown that the energy levels, correct to the second order, are the 
eigenvalues of an observable involving only the spin variables, namely 


Ey + 28(S3— NASH, A SS) = Oe eee (2) 


In these expressions f is the Bohr magneton, H the external magnetic field, 
\ the spin-orbit coupling constant and A,, is a real, symmetrical, positive definite: 
tensor defined in terms of the matrix elements of the orbital angular momentum 
L by 
(O|Ljla) (|L5|0) 
Ne od : Fi f df Alene 3 
je n#0 (ei) ( ) 


The tensor indices 7, 7 refer to Cartesian coordinates and the summation. 
convention is assumed for them. The indices 0,....7,.... refer to the orbital 
levels, 0 being the lowest; summation for these will be written explicitly. 
One can readily interpret each term in (2). The first, Eo, is the unperturbed. 
energy. ‘The second 
28S —AAy) SiH; 


is the magnetic energy of a spin system with a (in general anisotropic) g-factor- 
represented by the tensor 


SG = 20, NG) eee eee (5) 
in the external field H. ‘The third, 
—MAS)S;, a te etn ie. wee (6) 


1s the second-order contribution of the spin-orbit coupling. If A,; is isotropic: 
(i.e. a multiple of 5,,) it merely represents an equal downward shift of all the. 
2S +1 levels; if not, there is in addition a splitting of the levels, even in the 
absence of an external magnetic field. The last term 


7A EGE pity wtsiny ee eae ee (7) 


which is spin-independent and quadratic in H, corresponds to a temperature-. 
independent paramagnetic susceptibility. 
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§2. THE GENERAL CASE 

The perturbation procedure will first be described in terms of a more general 
system. For purely descriptive purposes the terms ‘orbital levels’ and ‘spin 
multiplicity’ will be retained, though the argument could equally be applied to 
systems whose unperturbed levels were characterized by other than orbital 
variables and in which the degeneracy arose in other ways. The unperturbed 
Hamiltonian will be denoted by “°, with eigenvalues Ey,...., E,,....  Eigen- 
states of #°,will be denoted by |0),...., |7),.... Each set {|)} is a linear 
manifold of dimensionality equal to the spin multiplicity. The total Hamiltonian 
will be written 


TE TE ate Hoge ae WET id a ek ate (8) 
The eigenvalue equation is 
epi) maT song) aa atin agence (9) 
One now expands | ) in terms of eigenstates of #°: 
pee Dees fe) ee ee ae. ae (10) 
n#0 


For a solution of (9) whose eigenvalue is approximately Z, |) is of zero order, 
while the |) are of first order. Part of the problem is to determine the appro-- 
priate |0) (secular equation). The eigenvalue equation (9) can be rewritten 


(OF +99] 04 DB E+) |(ny=00-- 20... (11) 
n#V 


The left hand side of (11) will now be decomposed into eigenstates of #°. 
This is conveniently done by means of the projection operators P,, of the linear’ 


manifold {|7)}. Operating on (11) with Py, P,,, one has 
ipa POS Portlny=0. » ates. (12) 
n#0 
(E, -E)|n) +P,3[0)+ EP, |my=0.  ..... (13) 


m0 
Equation (13) shows that |”) is indeed of first order in #1, and allows one to 
solve for it in first order as follows: 


oe 
=——* Se ae SEC Satire 14 
Jn)= = 10) (14 
Substituting (14) into (12) one has, to second order, 
Pot P gt s 
Het = a (ae Paes acer 15 
(2 + Py ep (FB) E)|0) 0, (15) 
or, in a more symmetrical form, 
y /1P 1p 
(z, D2 te i eal ~E)|0)=0. ee (16) 
nF (E,, — Eo) 


E, has been written for E in the denominators of (15) and (16) since the difference: 
involved is of third order. (16) is the required second order secular equation. 
determining |0). It shows that to this order the eigenvalues F are the eigen— 
values of the operator 


~ PWR, HP, : 
a opp eeate Ser Ate 8 eel 
H =E,+P,H'P, Pee aNes <3 (17) 


which.is Hermitian and operates entirely within the manifold {|0)}. 
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One may also determine the eigenstate to the second order. A straight- 
forward substitution of (14) into (13) gives |”) correct to second order, from which 
it is possible to give the following expression for | ), correct to second order, 


Deval nemee eo ees (18) 
where . : 
r=( eS Pa oe Se PH Pinte 
nZ~O (E,, — Eo) mA#z0n40 (E,, — Eo)(Em — £o) 


so; Pt Py) Ps 
nZ0 (E,,— Eo)? 0 


‘To derive this, 1/(E,,—E) has been expanded as 1/(E,, — Ey) +(£—£)/(£,, — 2)? 
.and (15) has been used to eliminate Z. The last term in 7, with the factor 3, 
has been included to ensure the normalization 


(\Y=40(0) ie: = Pes eee (20) 


‘The factor P, on the right of the definition of T is not essential, but as it restricts 
T to the manifold {| 0)} it has some advantages. 
As a corollary it may be seen that the expectation value of any observable € in 
the state |) can be written 


(A= Ole [05, = = eee (21) 


sie Pe ee A (22) 
A similar result is also true for the matrix elements connecting any two states 
|) belonging to the unperturbed level E,. In calculations dealing with the 
system in its lowest orbital level one may therefore use € and the states |0), 
‘instead of € and the states |); this frequently has practical advantages. One 


-example is the operator H, which is related to # in this way. Another, which 
will be discussed below, is the magnetic moment, — 6(L + 2S), in the paramagnetic 
problem, which can be expressed in terms of spin variables alone, the orbital 
moment being linked to the spin by the spin-orbit coupling. 


where 


$3. THE PARAMAGNETIC PROBLEM 

These general results may now be applied to the paramagnetic problem, for 
which the perturbation Hamiltonian has the form (1). The Schroedinger 
equation for the unperturbed problem does not contain spin dependent terms, 
but only the kinetic and electrostatic potential energies, both of which are real 
operators. ‘I'he wave functions for the states |”) may therefore be chosen to be 
real. From this it follows that the matrix elements of L, which is represented in 
the Schroedinger representation by the operator —7=rad/dr, are purely imaginary, 
and, since L is Hermitian, diagonal elements vanish and off-diagonal elements 


-are antisymmetrical ({m|L|]n)=—<n|L|m)). One therefore has 
P) #1 Po= P,[AL.S + B(L + 2S).H]P, 
= 2BS.H Pi ok ee eee (23) 


Py HP, HP, = (0| L;|n)<n| L;|0)(AS; + BH,)(AS,+BH,)P,.  (n40). 
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Hence 
a , OL; 1; 
H =, +28. H— y UE ENO Oras. + pH, 


n#0 (E,, — £,) 


(The factor Py has been omitted since it is to be understood that operates. 
within {|0)}.) On simplification and the introduction of the tensor A,; defined 
by (3) this reduces to (2). That A,; is real and symmetrical follows from the fact 
that (0| L|) is imaginary and equal to — {n| L|0). ‘That it is positive definite, 
provided og is the lowest eigenvalue, is also easily proved. 

One may also derive expressions for the effective magnetic moment. In this 
connection one seldom requires the second-order correction, which is small 
(provided A is small compared with the E,,—E,, which is the condition for the 
perturbation procedure to be a good approximation in this problem); it is therefore 
sufficient to quote the first-order results. They are 


FL SS) AROS ae) 2 yea a (26) 
ie Yad tet ov, ie Be signe ee (27) 
DOS) one ANS ORAM cote: (28) 


They exhibit the temperature independent susceptibility and the contribution of 
the spin-linked orbital moment to the g-tensor. 
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ABSTRACT. ‘The paramagnetic resonance from the nickel ion in nickel fluosilicate has 
been investigated at a number of frequencies and temperatures. Jt has previously been 
suggested that the ion is in an electric field of trigonal symmetry. ‘The results obtained are 
in good agreement with this hypothesis. It is also found that the magnitude of the splitting 
of the levels varies with temperature. 


Siz INTRODUCTION ANDI CRYSTAL VUOGRAPEHY 


URING the last two years a new means has been developed for elucidating 

the properties of paramagnetic ions in crystals. This is. the method of 

paramagnetic resonance, which, in principle, is very simple. A 
magnetic field is applied to the paramagnetic ions, so that the energy levels undergo 
a Zeeman effect. The ions are also subjected to a high-frequency magnetic field 
at right-angles to the steady field, so that transitions may be induced between the 
Zeeman levels when they have the appropriate energy separations, provided that 
the transitions are allowed. The consequent absorption of energy therefore 
shows a series of maxima as the static magnetic field is varied. By repeating the 
experiment for a series of different values of the frequency of the high-frequency 
field, it is possible to plot out the energy levels as functions of the static field.. 
- In particular, one can estimate the splitting of the energy levels in the absence of 
magnetic fields. 
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The situation is frequently complicated by the presence of several non-equivalent 
ionsinthe unit cell; here each type of ion contributes its own absorption spectrum. 
The fluosilicates of the iron group salts constitute a particularly simple case, 
however, since the unit cell contains only one paramagnetic ion. These salts 
form hexagonal crystals whose structure is known from the work of Pauling (1930) ; 
the paramagnetic ion is surrounded by six water molecules which form a slightly 
distorted octahedron. These give rise to an electrostatic field acting on the ion 
which is of nearly cubic symmetry but on which, owing to the distortion, a small 
trigonal component is superimposed. The distortion is in the direction of the 
hexagonal axis of the crystal, so that the axis of the trigonal field coincides with 
this crystalline axis. 

The behaviour of the Ni+* ion [(3d)§, 3F ; A= —335 cm] under the influence 
of crystalline fields was investigated theoretically by Schlapp and Penney (1932). 
A cubic field splits the F state into two triplets and a singlet; the singlet lies 
lowest and the triplets are of the order of 10,000cm™ higher in energy. The 
lowest level has a further threefold degeneracy because of the spin. When the 
combined effects of the non-cubic part of the crystalline field and the spin-orbit 
coupling are included, the spin sublevels undergo a small splitting of the order 
of lcm. At room temperature and below, only the spin triplet of the lowest 
orbital level is appreciably populated. The higher orbital levels are important 
as they enter into the perturbation calculations. 

Some of the magnetic properties of nickel fluosilicate have already been 
investigated by Becquerel and Van den Handel (1939), who measured the 
temperature variation of the Verdet constant. ‘Their experimental results were 
interpreted by Opechowski (Becquerel and Opechowski 1939), who investigated 
_ theoretically the behaviour of the Nit++ ion under the influence of a cubic field 
with a trigonal distortion. Opechowski found that in the absence of a magnetic 
field the basic spin triplet is resolved into a doublet and a singlet. The experi- 
mental results indicated that the singlet is 0-30 cm! above the doublet, and that 
the levels diverge with a g-value of 2-25, when a magnetic field is applied parallel 
to the hexagonal axis. By making certain assumptions about the form of the 
crystalline field, Opechowski obtained some information about the position of 
‘the higher orbital levels. 

The paramagnetic resonance experiments described in this paper cover the 
‘temperature range 14°-195°K. Some necessary extensions of the theory of 
Opechowski have been made by one of us (K.W.H.S.). The experimental 
results support the theory, and they reveal an interesting temperature variation 
-of the position of the energy levels. 


§2. EXPERIMENTAL METHOD 


‘The experimental technique will be described here only in outline, as it has 
already been discussed elsewhere (Bagguley, Bleaney, Griffiths, Penrose and 
Plumpton 1948). A single crystal of the salt under investigation lay at the bottom 
of a cavity resonator, in which oscillations of fixed frequency were excited by 
means of a klystron oscillator. A variable external magnetic field was applied 
at right angles to the radio-frequency magnetic field to which the crystal was 
subjected. The absorption of radio-frequency energy by the specimen reduced 
the Q-value of the cavity and was registered as a diminution in the power transmitted 
through it, which is proportional to QO. 
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At wavelengths of 3cm. and 1-2cm. the resonator was a cylindrical cavity 
excited in the H, mode. For longer wavelengths, quarter-wavelength concentric 
lines were employed. The resonator was immersed in a suitable constant- 
temperature bath, e.g. solid carbon dioxide in acetone, or liquid oxygen. 

It was necessary to use crystals of different sizes for the different wavelengths 
in order to obtain a suitable intensity of absorption—the 10cm. crystal had a 
volume of 0-5 cm*, while the one for 1-2cm. occupied 0:05 cm?. In some cases 
it was possible to use two crystals under the same conditions to afford a check. 


§3. EXPERIMENTAL RESULTS AND DISCUSSION 


In most of the experiments the magnetic field was applied parallel to the optic 
axis (c-axis and hexagonal axis) of the crystal, because the magnetic behaviour is 
then particularly simple, the energy levels diverging linearly (see § 4) 


W,=+2,8H—6; W,=—g,BH—5; LEAVE 


where 6 is the initial separation, g, the ‘g-value’, and B is the Bohr magneton. 
It is assumed that the doublet is lower than the singlet, so that 6 is positive, in 
agreement with the findings of Becquerel, Van den Handel and Opechowski. The 
paramagnetic resonance experiments give no evidence on this pcint; they are 
consistent with either the doublet or the singlet being the lower. 

It can be shown that in the present case transitions are not allowed between 
W,and W,. The transitions W,<—- W, and W,<— W, are permitted, however, 
and they give rise to absorption peaks for values of the magnetic field given by 


|¢, BH —3|=hAv and |g, BH +8|=hv. 


The absorption spectrum was investigated at wavelengths near 1-25, 3-2, 7 and 
9cm. Theresults foratemperature of 90° kK. are presented graphically in Figure 1. 
They confirm the expected linear divergence of the levels, and yield the value 
0-17 cm“ for the initial splitting. The mean g-value is 2-25,. 

When the paramagnetic resonance experiments were repeated at other 
temperatures, it was found that the initial splittings changed considerably. 
The values for g and 6 for different temperatures are given in the Table. 


Temperature (° K.) 195 90 60 20 14 
Sine 2:29 2:26 — 2-29 Eee, 
6 (cm?) 0-32 0-17 0-14 0-12 0-12 


Evidently g is independent of the temperature, within the limits of experimental 
error; 5, on the other hand, varies by a factor of 2-5 in the range of temperature 


investigated. . 
In attempting to explain these phenomena one must anticipate some of the 
conclusions reached in the theoretical section, where Opechowski’s calculation 


is extended. He has shown that 

De (Aper tate aah Mele) FES ae (1) 
and 
| Eu a 2(1 —AA), Ska tsMauie (2) 


here A and B are parameters which depend in a complicated way on the electric 
field acting onthe ion. The spin-orbit coupling coefficient, A, is known from the 
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work of Laporte (1928) to be —335cm™. Combining this with our value for 


5 at 90° (0-17 cm*) one finds 
A=(2~—g,)/2A =(2-26 —2)/2 x 3353-9. 10°-*cm., 


while 


(A —B) =8/2 =0-17/(335)2 1-5. 10-8em., 


which is about 250 times smaller than A. 6 thus represents the difference 
between two much larger quantities which are determined by the field parameters. 

- If the electrostatic field were purely of cubic symmetry, 6 would be Zet0; 
this, however, is not so, and the actual value of 6 depends in a rather complicated 
way on the magnitudes of the cubic and the trigonal terms in the electric field 
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Figure 1. A paralle! to hexagonal axis. 
Temperature 90° k. 


acting on the ion. It is not determined solely by the magnitude of the trigonal 
field. ‘The variation in § with temperature could thus be due to changes in the 
internal field due to thermal expansion of the crystal. 

In another series of experiments the magnetic field was applied at right angles 
to the hexagonal axis, to provide a further test of the theory. The energy levels 
do not now diverge linearly with magnetic field but are given by 


Wi=—8;  Wa= —3B— (844g 6H]; Wy= — Hb + ("+ 49,262), 


where g, =2(1—AB). As we have already mentioned, (A — B)/A~1/250, so that 
(£1 —-£,)/2u~1/250. Thus for practical purposes we can say that g,=g,. 
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In strong magnetic fields (¢4¢8H>5) the energy levels can be written as, 
approximately, 


W,=-5; W,= —36+ 8, 8H; W,= —46—g, PH. 
The selection rules for transitions between these levels induced by an oscillating 
magnetic field at right angles to the steady field are that the transitions W, <— W, 
and W, <> W, are allowed, but that the transition W,<— W, will be very weak. 
Thus at a given frequency v two strong absorption lines will be observed at values. 
of H such that hv=g,BH— 36 and hv=g,BH +38. 

Here again it was possible to plot the behaviour of the energy levels (Figure 2). 
The curves which have been fitted were calculated using for g, the value found for 
g, in the previous experiments. The curved part could not be properly verified 
on account of the low intensity of absorption at the long wavelengths concerned. 
The available points, however, lie satisfactorily on the curves. 


§4. THEORY 

In order to explain the results of Becquerel and Van den Handel, Opechowski 
considered the energy levels of the nickel ion in the presence of asteady magnetic 
field directed along the hexagonal axis. ‘To explain the results given graphically 
in Figures 1 and 2 it was thought desirable to extend his work and to calculate 


W, 
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Figure 2. H perpendicular to hexagonal axis. 
Temperature 90° k. 


the energy levels as functions of the external field and also of its orientation with 
respect to the hexagonal axis. This requires a somewhat different treatment 
from that given by Opechowski, as the energy levels do not, in general, diverge 
linearly with increase of magnetic field, and also the changes in the energies are 
often large compared with the separation of the levels in the absence of the field... 
PROC. PHYS.-SOC. LXIII, I—A 3 
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The effect of the electric field on the orbital states is the same as that given 
by Opechowski, and his work is applicable thus far. To find the combined 
effect of the spin-orbit coupling and the magnetic field, for any direction, it is 
necessary to take them together and to treat them as a perturbation on the lowest 
spin multiplet. A straightforward application of first order perturbation theory 
shows that, while there is a contribution from the magnetic energy, the contribution 
from the spin-orbit coupling is zero. In proceeding to the second order change in 
energy it would not, in general, be correct to use the zero order states which make the 
magnetic energy diagonal, except when the change in energy due to the magnetic 
field is large compared with that due to spin-orbit coupling. When the magnetic 
field is along the hexagonal axis, it is found that the zero order states which make 
the magnetic energy diagonal are also those for which a straightforward application 
of perturbation theory is valid. | As Opechowski was only interested in this case, 
he avoided the difficulty given above. 

A method suitable for dealing with this type of problem has been given by 
Van Vleck and is described in a paper by Jordahl (1934). Another method is 
described in a note by Pryce (1950). 

The perturbation is taken as 


AL.S— BH. (L +28), 


where A is the spin-orbit coupling coefficient, WL is the orbital angular momentum 
operator, 4S is the spin angular momentum operator, f is the Bohr magneton, 
and H is the steady magnetic field. Using the second of the above methods, 
to second order, it is found that the energy levels of the lowest spin multiplet are 
the same as the eigenvalues of 


—2BH.S— A(02S.2—2\BH_S. + 62H,2) 
— BIA(2= S,") — 2B HS, +, 5,) p(y ea) 


where 4 and B are constants which, in Opechowski’s notation, are given by 


ee ae i ae eee 
Rima 1S Pe aes 


and where the hexagonal axis is taken as Oz. The spin-independent term 
— AB*H,? — B20? + BH? + H,?)], 


can be dropped, since it corresponds to a displacement of all the levels, though 
it is of interest in that it shows that there is a displacement proportional to the 
square of the magnetic field, which will give rise to a constant term in the magnetic 
susceptibility. By writing 2(1—AA)=g,, 2(1—AB)=g, and (A—B)A2=8, the 
energy levels are the eigenvalues of 


= Ose = PLL gS; c& Bg (HS, ad S, Js 


The significance of 6 is apparent since it represents the separation of the levels 
in zero magnetic field. 


Y 


In order to find the selection rules governing transitions between the 
components of the lowest spin multiplet, it is necessary to determine their states. 
A good approximation to the actual states will be given by taking them as 


products of the lowest orbital state with the eigenstates of the operator (3) 
given above. 
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‘The states on which measurements are made have been obtained for the two 
directions 6=( and @=7/2, where 6 is the angle between the direction of the 
magnetic field and the hexagonal axis. It is convenient to take the hexagonal 
axis as axis of quantization, so that S, is diagonal with eigenstates |1>, |0> and 
|—1). For the direction @=0, the energy levels and states are 


W,=g,PH—56; Peas 
Wee yell 0) ly, 
Wo=0 |0), 


where the orbital state has been omitted as it is the same for each level. The 
matrix elements of all components of L are zero between these states, so the 
selection rules for an oscillating magnetic field perpendicular to Oz are such that 
the transitions W,<— W, and W,<—> W, only are allowed. 

For the direction 6=7/2, (Ox), the energy levels and states are 


fotent (I1)=1=19} 

Wee VO ag, PP), 2,RH1) +) —1) pa 2|0); 

W,=— 319+ V(8'+4e°R°H*)];  g, BH{|1) +] —1)} + V2W,|0). 
It is to be noted that the states are not normalized. The oscillating field may now 
be along either Oz or Oy, but in either case it is found that the allowed transitions 
are W,<— W, and W,<— W,. 


For a general direction of the magnetic field the energy levels are given by the 
roots of the secular equation: 


i : 
—d5—fe,Hcosé@-W - <9 Bg. H sin 6 0 
V 
! / : ee aie 
- Borns sin? —W _ 72 PB. sin =f) 
0 — 5 AeyH sin 6 —6+fe,Hcosé—-W 


‘This is a cubic and does not in general factorize. 


§5; CONCLUSION 

The results of the paramagnetic resonance experiments are entirely consistent 
with the crystalline field theory developed by Opechowski and extended in this 
paper. The magnetic properties of nickel fluosilicate have also been investigated 
by two other methods. Becquerel and Van den Handel (1939) and Becquerel and 
Opechowski (1939) measured the temperature variation of the Verdet constant 
in the temperature range 1-5°k.-290°K., and found that the magnetic moment of 
the ion is 2:25 Bohr magnetons when the magnetic field is parallel to the hexagonal 
axis of the crystal. The deviations of the Verdet constant from a 1/T law enabled 
them to compute a value for 5. ‘Their measurements are sensitive to ) only at 
the lowest temperature; hence their value (0-30 cm) refers to the helium region. 
On the other hand Benzie and Cooke (1950) find 6=0-15—0-16cm ? at similar 
temperatures. This can be reconciled with our own measurements at higher 


temperatures. 
2 
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The experimental work described in this paper was done in the summer of 
1948 by Dr. Penrose, working at the Clarendon Laboratory, Oxtordae ie 
subsequently moved to Leiden where he intended to continue these experiments. 
at helium temperatures. Owing to his illness, from which he has recently died, 
the experiments were not done, but it has seemed worth while to publish these: 
results, though they are not as complete as originally intended. 


Note added in proof. Holden, Kittel and Yager (Phys. Rev., 1949, 75, 1443): 
have reported a measurement at room temperature and 1:25 cm. wavelength on 
NiSiF,, 6H,O. The data for magnetic field parallel and perpendicular to the 
trigonal axis are not quite consistent with the assumption of trigonal symmetry, 
but yield the values 8 =0-49—0-52 cm~', g=2.29—2:36. ‘The value of 0 falls 
on a smooth curve with our measurements at lower temperatures. 
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The Theory of Magnetic Resonance-Line Widths in Crystals. 
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ABSTRACT. A theory is developed for calculating the mean displacements and mean 
square widths of resonant absorption lines in crystals for which the spin-lattice relaxation 
time is long compared with the spin-spin relaxation time. The theory is applied to a 
number of different cases ; in particular, the effects of nuclear hyperfine structures and 
exchange forces in ionic resonances are discussed at length. It is shown that the theory 
holds, provided that the temperature is high compared with the Curie temperature. : 


Si LN ERO DWEeDION 


WE recent development of the method of magnetic resonance has provided 

a new and powerful method for studying the properties of magnetic systems 

in solids and liquids. Its application to crystals containing paramagnetic 
ions has led to a considerable increase in our knowledge of the fields acting on the 
ions, most of which had previously been obtained from magnetic susceptibility 
measurements. 

The principle of the method is very simple. If an ion possessing a magnetic 
moment, and therefore having a degenerate ground state, is placed in a steady 
magnetic field, the degeneracy is lifted and the levels undergo a Zeeman splitting. 
An oscillating magnetic field will be able to induce transitions between the Zeeman 
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levels if it has the right frequency and if the transition is not forbidden. Thus, 
if the frequency of the oscillating field is varied, the absorption of energy from it 
will show a series of maxima as the frequency passes through the characteristic 
frequencies of the ion. For a free ion there will be only one characteristic 
frequency, but when the ion is not free there may be more than one frequency of 
absorption. 

For paramagnetic ions in fields of order several thousand gauss the characteristic 
frequencies usually fall in the microwave region, while for magnetic nuclei in 
similar fields they fall in the radio-frequency region. The experimental methods 
have been described by a number of observers (Bagguley et ai. 1948, Bloembergen, 
Purcell and Pound 1948) and will not be given here. 

As is well known (Van Vleck 1932), ions in crystals are usually in strong electric 
fields. For a paramagnetic ion the fields distort the orbital motion of the electrons 
and are also coupled to the resultant spin through the spin-orbit coupling. The 
effect of this is that the ground state of the ion in a crystal is often very different 
from that of a free ion and its magnetic properties are correspondingly altered. 
Usually the orbital motion 1s ‘quenched’ and the magnetic moment of the ion is 
anisotropic with a value near that of the ‘spin-only’ value. 

Magnetic resonance experiments have confirmed these ideas and have given 
accurate values for the magnetic moments of ions in different crystals. They have 
also shown that the widths of the absorption lines (usually of order several hundred 
gauss) are much greater than the natural widths which would be expected if the 
ions could be regarded as independent systems placed in suitable steady eiectric 
and magnetic fields. In other words, the models which have usually been used in 
calculations of the positions of the resonance lines and the magnetic susceptibilities 
are not good enough approximations to give the line widths. 


$2) A GENERAL SURVEY OF THE PROBLEM 


In a theory of Jine widths in crystals it is necessary to consider the interactions 
between the ions themselves and also between the ions and the thermal motion of 
the lattice. A general treatment on these lines has not yet been given. In this 
paper an account will be given of a method which can be used when the interaction 
between the ions and the lattice is very small compared with interactions between 
the ions themselves. Under these conditions a good approximation is given by 
regarding the magnetic ions as located at fixed positions and interacting through 
their couplings with each other. The thermal motion will modulate the inter- 
actions and the approximation is equivalent to assuming that the extra width from 
the modulation is small compared with that which would be found were there no 
thermal motions. This is not an unrealistic approximation, as there are many 
crystals in which the widths of the resonance lines are determined almost entirely 
by the static interactions. These are crystals in which the spin-lattice relaxation 
time is long compared with the spin-spin relaxation time. 

Before giving the quantum mechanical treatment of the line widths under these 
conditions, it is instructive to consider the way in which broadening arises, froma 
semi-classical point of view. Each ion is regarded as a gyroscopic magnet; 
that is, its angular and magnetic moments are parallel. We consider a system of 
such magnets, each located at a fixed point in space and all in a uniform steady 
external magnetic field, directed along Oz. Each magnet will precess about Oz 
and can be regarded as equivalent to a magnet fixed in the direction of Oz together 
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with a magnet rotating in a plane perpendicular to Oz. Ifa magnet Ais considered, 
the steady field on it will not be equal to the external field but will be the resultant 
of the external field and an internal steady field from the fixed magnets associated. 
with all the other magnets. A will therefore precess about this resultant field with 
a frequency equal to yH, where y is its gyromagnetic ratio and H is the resultant 
steady field. The value and orientation of H will vary from one magnet to another, 
with the result that there will be a spread in the precessional frequencies. ‘This 
will appear as a broadening of the resonant absorption line. It 1s to be noted that 
this type of broadening does not depend on the magnets all having the same 
gyromagnetic ratios, and may be termed ‘steady field broadening’. 

It is next necessary to consider whether the fields set up by the rotating 
magnets can have any broadening effects. If two magnets have the same precess- 
ional frequencies, the rotating field from one will be at the right frequency to induce 
transitions in the other (the transition frequency is equal to the precessional 
frequency), whereas if the two magnets have different precessional frequencies, the 
rotating fields will have little effect. This shows that magnets with the same 
resonance frequencies will tend to induce transitions in each other and reduce 
their life-times in given states. This will appear as a broadening of the resonance 
line. There is no way in which the total energy of the system can change, so these 
transitions will only occur if they are energetically allowed. They can be pictured 
as transitions in which one magnet gives energy to another. This type of 
broadening will be referred to as ‘resonance broadening’. 

In the quantum mechanical treatment of the line widths the expressions. 
obtained are usually fairly complicated, but it is often found chat the matrix 
elements which occur can be interpreted as arising either from resonance 
broadening or from steady field broadening. Further reference to this will be 
made after the quantum mechanical treatment. 

The semi-classical treatment has been given under the assumption that the 
only type of coupling between the ions is magnetic dipole-dipole coupling; in 
actual crystals there may also be exchange couplings. There is no difficulty in 
taking exchange into account in the theory to be given; the interaction used can 
include both exchange and magnetic coupling. 

The method used is a perturbation method in which the energy of interaction 
is taken as a perturbation to the energy of an uncoupled system of ions. For this, 
it is convenient to assume that the energy levels of a single ion, under the influence 
of the crystalline electric and the applied magnetic field, fall into a low-lying group 
with separation small compared with RT, and that the other levels are so far 
removed in energy that their populations are negligible at the temperatures used. 
These higher levels play no part in the theory and it is convenient to forget them 
and assume that the ions have only the low-lying levels. This restriction is not 
severe, as it usually happens that the levels are arranged in this way. 

It is found that, provided RT is large compared with the energy differences 
between the states in the low-lying group, the mean square width of the absorption 
lines is effectively independent of temperature, and that the centre of the broadened 
line is undisplaced from the position it would have were there no inter-ionic 
couplings. At temperatures for which RT is of the same order of magnitude as 
the energy separations, both the displacement and the line width may be tem- 
perature dependent, and there is also a dependence on the shape of the crystal, 
an effect closely connected with the demagnetization factor (see below). 
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The theory is also applicable to calculating the widths found in nuclear 
resonance experiments, and it is useful to have a notation which can refer to either 
ions or nuclei. ‘To this end we shall use the word ‘dipole’ to mean either a 
nuclear or an ionic magnetic moment. The magnetic coupling between dipoles 
will be referred to as ‘dipolar coupling’. 

The oscillating magnetic field is usually at right angles to the steady field in 
resonance experiments. It is not necessary to assume this in the theory and the 
formulae hold for any angle between the fields. 


§3. THE DERIVATION OF GENERAL FORMULAE FOR THE MOMENTS 


In this section, formulae of general application are obtained for the total area 
of an absorption line, and its first and second moments about the position it would 
have were there no inter-dipole couplings. From these quantities the mean 
displacement and mean square width can be found by taking the appropriate ratios. 
The absorption curve is taken as a curve in which the absorption A(v) at frequency v 
is plotted against vy. A(v) can be measured in arbitrary units, since in taking the 
ratios any common factors will disappear. The method is readily extended to 
calculating higher moments, but the relevant formulae will not be given, as they 
are, in principle, easily derived, though cumbersome. 

It is supposed that the crystal contains dipoles which have a number of possible 
transitions, one of which occurs at a frequency »v, and it is the width of the line at v 
which has to be calculated. In the absence of any coupling, the energy levels of 
the dipole system will consist of a large number of highly degenerate levels, and 
there will be many pairs of levels with separations hv. One such pair is chosen, 
and their energies are denoted by E, and Ey, so that hy=E,—E,. The manifold 
of degenerate states in E, is described by a projection operator P,, and that for 
E; by P;. If the effect of a coupling is now considered, the degeneracies are 
largely removed, and instead of the levels E and E, giving a sharp line at v, the 
states which have come from them will give transitions which are spread about v. 
This process will occur for all pairs of levels with initial separations Ay, and in 
determining the total width it is necessary to average over all such pairs, with 
appropriate weighting factors. It is assumed that all the width is due to this 
cause, and that there are no contributions from states which were initially at 
separations other than hy. If such contributions were appreciable the spectra 
observed would not consist of the distinct lines found, but, instead, the lines would 
be broadened into each other. 

The inter-dipole coupling is represented by a term W in the Hamiltonian of 
the system, and the zero order states in the perturbation calculation are chosen to 
be those eigenstates of the total energy which make W diagonal in the manifolds 
P, and P;. Then, if |a@) represents such a state in P, and |b) a state in Ps, 
the change in energy of |a) when W is introduced is, to a first approximation, 
{a| W|a), and that of |b) is (b] W|d). 

In order to study the resonant absorption when an oscillating magnetic field of 
frequency near to v is applied to the crystal, it is necessary to consider its contri- 
bution to the Hamiltonian of the system. This will be denoted by M, and it can 
be regarded as inducing transitions between the various states of the dipole system. 
It will cause transitions between states |2) and |b) if the frequency of oscillation 
is v+{(b| W|b)—<a| W|a)}/h, with a transition probability proportional to 
| <a|M|b)|?. The probability that the crystal is in the state |a) is different from 
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that for state [6), and it is this difference which makes the rate of absorption of 
energy greater than the rate of stimulated emission. ‘When allowance is made for 
this, the spectrum from states |a) and |), apart from constant multiplying factors, 
can be expressed as: 


(exp (—E,/RT) —exp(—B,/RT)}(Ey—E,)| <a] M| 6)? 
x 5{v' —v—{<b] Wb) — ¢a| Wa)}/h}/Z, exp (— ERT), 


where E, is the energy of the state |a), EZ, is the energy of the state |b), Tis the 
absolute temperature, Rk is Boltzmann’s constant, the summation is to be taken over 
all possible values of E,, v’ is the varying frequency, and 6(v’) is a function of v’ 
which is sharply peaked at v’ =0. (v’) is introduced to represent the shape of the 
absorption line. If this expression could be summed over all pairs of states such 
as|a) and|6), it would give the shape of the complete absorption curve as a function 
of the frequency v’. This, however, is not possible, but by making certain approxi- 
mations it is possible to obtain a good deal of information about the shape of the 
curve, for it appears that the area, mean centre, mean square width, and higher 
moments can be written in invariant forms which can be evaluated without 
knowledge of the zero-order states. It is not practicable to calculate more than 
the first few moments, as the labour involved increases rapidly with increasing 
order. 
The expression 


{exp (—E,/RT) —exp(—£,/RT)}(E,—£,) 
can be written as: 
exp(—E,|kT) exp (— (al W]a)/RT)[1 —exp {—(hy + (|b) — (al Va) RT} 
[hy + (6| W|b)— <a] W] a>], 
where E, is the energy of the manifold P,. As <«b]W|b)—<a|W|a) is small 


compared with hy for all the pairs of states which make appreciable contributions 
to the width, the error will be small if the above expression is replaced by 


exp (—E,/RT)exp(— (a|W|a)/RT)hv{1 —exp(—hv/RT)}. © 


This approximation involves no assumption about the magnitude of the 
temperature. If it is now assumed that the temperature is much greater than the 
changes in energy caused by the interaction, a further approximation is to replace 
exp {—(a|W|a)/RT} by unity. The error introduced by this approximation 
will be considered below. The absorption curve for the states |a@) and |b) can 
then be written as: 


exp (—&,/RT)hv(1 — exp {— hy/RT})| (a| M|6)/? 
x o[v' —v —{¢b| Wb) — a] W| a)}/h] /X, exp (—E,/RT). 


In finding the mean displacement and mean square width of the absorption line, 
it will ultimately be necessary to consider the ratios of expressions derived from 
the above formula, so it is convenient to drop the factor 


hyv{1 — exp (— hv/RT)}/X, exp (—E,/RT), 


as it does not depend on the particular choice of the states |a@) and |b). The total 
area of the absorption curve is then found by summing the areas given by all pairs 
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of states such as|a) and[b). If the integral of 8(v’) for all values of v’ is taken to 
‘equal unity, and if the summation is first taken over all pairs of states coming from 
manifolds Ay apart, the area can be written as: 


diy, p EXP ( HIRE) Spur P,MP,M, aesteme (1) 


where the summation is to be taken over all values of « and 8 such that E,—£,= hv. 
In the same way, the first moment about the undisplaced position is found to equal 


2, exp (—E,/RT) Spur (P;>WP,MP_M — P, WVIAC RAY 8 eae cai: (2) 
-and the second moment: 


=,,¢ exp (— E,/RT) Spur (P;WP,WP,MP,M —2P,WP,MP,WP,M 
SE WRWE MPM 9 es (3) 


The ratio of (2) to (1) gives the mean displacement of the line from Av, and the 
ratio of (3) to (1) gives the mean square width. In taking these ratios, it is seen 
that the various constants which have been dropped would disappear had they 
.still been present, a result which justifies their earlier removal. 


04 ONE yeh Or eDTP OI 


The above formulae are most easily applied to the case of a crystal in which 
there is just one type of dipole present, and in which the energy levels of the single 
-dipole are such that there are only two states with energy difference hv. It is 
supposed that each dipole has R states which are denoted by |a@,), |a.),....|ap), 
with corresponding energies @,, a, .... Gp, and that the transition of interest is that 
between |a,) and |a,), so that a,—a,=hv. The energies of the manifolds P,, can 
then be specified by the number of dipoles 1n the various states, and if N is the 
total number of dipoles in the crystal, a particular pair of levels E, and E, may be 


specified by the sets of numbers (1, 19, 13,....mg) and (m,—1,n,+1,n3,....7,) 
where 

EB = 1G, + eGg + Nyda + .... +Npap 
and 

Ey =(m,—1)a,+(m,+1)a,+ngag+.... +npap 


with Li=#n,=N. This replaces the summations over the «’s and fs by 
summations over the n’s. 

Before carrying out the summations, it is convenient to decompose M and W 
into their component parts, and to write M=X,m, and W=%,,_,W,,,, where m, is 
the coupling of the radio-frequency field to the sth dipole and W,,, is the inter- 
dipole coupling between dipoles p and g. In order to simplify the notation, it is 
useful to write expressions of the form 


CO ste NOY iG, as an 
as 


Cis eel Dad Adare) 


where the a’s are omitted altogether. This leads to no confusion as long as there is 
-only one type of dipole present, and it so happens that even when there are several 
types present, the final results can be written in this form without ambiguity. For 
dipolar and exchange couplings (which are of most interest), expressions of the 
form »,{,s|W,;|,s), being of the nature of spurs over the states of the jth ion, 
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are zero. It will therefore be assumed that all such expressions are zero. The 
summations can be carried out without this assumption, but it is convenient to 
make it, as it leads to a good deal of simplification. 

The results of the summations are that (1) and (3) give constant terms together 
with temperature dependent terms which are completely negligible if RT is large 
compared with the a’s, whereas (2) gives terms of order a,/RT. This means that 
for the temperatures normally used, the mean displacement of the line is small, 
and that the area and width may be calculated as if the temperature is infinite. 

It is readily seen that 

Spur P,m,P,m,=0 


for Spur Pm, Pott) = 255 o . . \G0, Cre Eat belt Ga eens 


and each term is certainly zero unless the ¢th dipole is in the state |1); but in this. _ 
case the operator Pm, changes it into the orthogonal state |2), and as there are no 
other operators acting on the ¢th dipole, these terms also vanish. It also follows. 
by symmetry that 

Spur P,m,P,m,= Spur Pym, P3m, 


and 
Spur Pit Patty = Lig ye... (ly Gs Oss | atte gt, | WoO roe) 
=e pe Al ml2) Gee oes | le eaee, 
n : (N—-1)! 
Kee (n,—1)! n,!75!....mg!’ 
so that: 


©,,9 exp (—~H,|kT) Spur P,MP,M=N| (1|m|2)Pexp(—a,/RT)A™, 
where 
A = d= exp (—a,/RT)=R. 


Coming now to the first and second moments, it is probably simplest to evaluate: 
(2) and (3) by writing them as summations of the form: 


xj, i hice, ODUL (P, W,;P,m,,P,m,—P,W,,P. m;,P3m) 


a 


and Xi 5,n,.... Spur (PW P,WiyP ym, Pam, —2P,W,;Pm,PsW Pam, 


CALNE) yy” ap 


+ P3W,;P,W).P3m,P,mq); 


and then to evaluate each term separately. This is quite long, so only the results. 
will be given. With RT large compared with the a’s, (2) vanishes to first order in 
ak T. No terms of type W,,W,,, survive in (3), and terms of type W,; W,;, Survive 
with factors of order a,/RT, compared with those of type W,;W,,, and can thus be 
neglected. ‘The value of (3) is then given by the expression: 
[1,1] Wy] 1, 1) — (2, 1] Wyl2, 1) — (2, 1] Wl, 292 
scafm29 preg; | tI 21 Mal 12) ~ €2,2] W122) +42, 1] Weyl, 290 
4 SUSE | [{1, 5] Wis] 1,8) — <2, s| Wy|2, 5)? | 
| +141, s| Wels, 1) — <2, 5] Wiss, 2)? J 


s>2 


If the summation over all 7 andj of the terms in the square bracket is denoted by X,, 
the mean square width is seen to equal X/NR. 
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For the sake of uniformity with the expressions found for the areas and second 
moments in more complicated cases, it is useful to write the above expressions in 
the forms: 


N 
Area D | <1] m|2)/. ne (4) 
Second Moment 
(41, 5] Wj] 1,8) — <2, s| W,,|2,s)}<2|m] 1) 2 
By | +B, 42,4] PrsWipPasl 1, R)CR| |) 
ee | — 2), 42,1] PosWizPos| 1,2) <k| m1) 6) 
oat iaaw | ye ) OScee an vars ie) 
| <1, s] Wijls, 1)— 2,8] Weyl, 29} 2|om| 1) P 
+2.) + Uy 42, 1] P1WiPi5| Rk, 1) ¢R| m| L) 
— Xi, x42, 1] PosWizPos| Rk, 1) <R| m| 2) 


Here, D denotes the total number of states of the crystal and the dash on the second. 
summation over s means that the values 1 and 2 of s are to be excluded. This. 
avoids the repetition of terms which are already included in the first sums. The 
operators P,, and P,, are to be interpreted as the projection operators of the states. 
of two dipoles corresponding to energies a, +a, and a,+a,. The possible values. 
of / and k are then severely restricted, and the whole expression is exactly equal to 
that given earlier. 


SSe “ENVGy APNOPISS: (un IDIUEONEAs 


In this section it will be supposed that the crystal contains two types, A and_ 
B, of dipoles, and that no transition frequencies of the dipoles coincide. As before, 
the states of the dipoles of type A (Rin number) are denoted by|a,),|a,),...-|@p), 
and those of type B (S in number) by |6,), |4,),....]s), with corresponding 
eNergieS G1, do,....@pz; 51, by,....bs5. The transitions between states |a,) and 
|a,) are considered, and it is again supposed that the partial spurs 4, ¢,s| W;;| , 5) 
are zero, where dipole j may be of either type. 

With the same approximations as before, the first moment is zero, and the others. 


are: 
Area p28) <1|[m,|2)|?. 
Second Moment 
1 
DA] ma|2)PX 
1 9 
+D yell] ms [2 PE 2it%I (1, s|W,,[1, 5) — <2, 5] W5|2, 5/75, 


where D is equal to the total number of states of the crystal, N, is the number of 
dipoles of type A, X is equal to the same expression as before with dipoles z andj 
both of type A, and the second summation over 7 andj is to be made with dipoles 7 
of type A and dipoles 7 of type B. 

The general expressions (4) and (5) will include these results if Vy is used for 
-N, and the summations over s are given a wider interpretation, so that s stands for 
states of both kinds of dipoles, and dipole 7 is always taken to be of type A, whereas. 
dipole j can be of either type. When dipoles 7 and j are of the same type, there is 
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-a factor 1/R?, and when they are of different types the factor has to be altered to 

1/RS. Certain of the terms which then appear have no meaning, since a dipole 
.can only be of one type, and it is necessary to take all such terms as zero. The 
-extensions which have to be made if there are more than two types of dipoles 
present, and there are no coincident lines, are readily made. All that is necessary 
is to let s range over the states of all the dipoles, and to change the factor 1/R® in 
-an obvious manner. 


S65) «COUN Cl Dit Nae Se INES 


It frequently happens in the case of nuclear resonance, and sometimes in 
paramagnetic resonance, that a number of transition frequencies for a single dipole 
coincide. When this occurs, there are resonance effects between the states (as 
discussed earlier), and the above formulae have to be modified. 

In this section it will be supposed that there is only one type of dipole present 
-and, with notation as before, it will be assumed that a, —a,=a,— 43, so that the 
transitions a,<—>a, and a,<—> az coincide. The area of the line is found to be 
equal to the sum of the areas calculated as if the lines were distinct and the mean 
displacement is, as usual, zero. ‘The second moment is found to be given by the 
sum of two expressions such as (5), one for the transition a,<—> a,, and the other 
for the transition a@,<—>a3. ‘The summations over / and k will not now be quite 
so restricted, and it is to be noted that the width of the line will depend on the 
matrix elements of the magnetic moments of the dipoles. If more than two lines 
coincide, there must be a term of type (5) for each coincidence. 

A similar result is found when the coincidence is of the form given by 
Ay — @, = 3 — Ay, So that the transitions a,<—> a, and a,<>a, coincide. The area 
is again equal to the sum of the areas calculated as if the lines were distinct, the 
mean displacement is zero, and the second moment is given as a sum of expressions 
like (5), with one for each transition. 

If other types of dipoles are also present, the formula for the area is unchanged, 
-except that the value of D is altered, and in the second moment it is necessary to 


let s range over the states of the dissimilar dipoles as well as the similar dipoles, as 
before. 


§7. LINES FROM DIFFERENT DTYPES OF DIPOLES COINCIDING 


In crystals which have two non-equivalent dipoles in unit cell, there are usually 
orientations of the applied magnetic field for which the two types of dipoles have 
the same energy levels. When this case is examined, formulae (4) and (5) are 
still found to hold, each coincidence giving its own contribution. In the cases 
previously considered, the summations over / and k give zero for dissimilar dipoles, 
but in this case they give non-vanishing contributions. 

It thus appears that formulae (4) and (5) can, with appropriate interpretations 
of the summations, cover many of the cases actually met with experimentally. 
‘The areas are always directly additive, and the second moments are also additive, 
with the difference that the actual expressions to be added depend on the problem 
under consideration. As the changes to be made always take place in the summa- 
tions over / and &, it is reasonable to interpret these terms as the contributions 
-arising from the precessional resonances described earlier. 
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§8. DEGENERATE AND NEARLY DEGENERATE LEVELS 
The external magnetic field is usually sufficient to ensure that the states of the 
dipoles are non-degenerate, and this has been assumed in the preceding work. 
It is possible, however, that for some values of the applied field two energy levels. 
coincide. ‘This happens, for instance, in chrome alum for certain orientations 
and magnitudes of the field (Bagguley and Griffiths 1947). It is necessary to 
consider this case separately. Let it be assumed that only two levels coincide 
and that they have energy a,. The basic states for this level may be chosen to be- 
any two mutually orthogonal combinations of the states in a, say |a)) and |a;). 
The expressions for the intensity and the width will depend, inter alia, on 
{a,|m|a,) and {a,|m|a;). It is convenient to choose as basic states those linear 
combinations which make (a,|m|a,;)=0. The results must be, of course, 
independent of the choice of basic states, and can be written in general form when: 
the result for the restricted choice is known. 
The area of the absorption curve is 


D2 || m|2yP. 
and the second moment is 
= id |m| 2p 
j 2| <2, 1| W,;|3, 1) — (2, 2] W,]3,2)+ <1, 2] W,]3, 1, | | 
x | X+E EY +1 42,1] Wy 3, YP +12, 1) Wel, 3dP 
t + |<2, 2| W|3, 3>P + [¢2Z, 2[ Wy] 2, 3)? a 


X is the same as previously defined for the transition a,<-+a,. The extra terms. 
giving the effect of the degeneracy are not given by (5) and have therefore been 
written out in full. 

If a dipole has two energy levels, say a, and a, which are very near to each other, 
the transitions a,<-~a, and a,<—>a, give two lines with a small frequency 
separation, equal to |a;—a,|/h. When, however, asystem containing such dipoles 
is considered, the lines are broadened by the inter-dipole couplings and it may 
happen that the widths are comparable with the separation | a3 —@,|/h, so that the 
lines overlap and can no longer be regarded as distinct. It is then incorrect to 
calculate the second moment as if the lines are distinct, for if this is permissible,. 
the result so obtained should be independent of |a,—a,|/h, and as this tends to. 
zero should go into (6). This does not happen. Similarly, it is incorrect to 
calculate the second moment as if the levels are degenerate, using (6), for if W is 
put equal to zero, the result is certainly wrong. In order, then, to deal with nearly 
degenerate levels, the theory has to be extended, as follows. 

Consider, first, the Hamiltonian of a single uncoupled dipole, and let its 
eigenvalues be aj, @, a3, dy, .... etc., where a, and a, are nearly equal. If the 
corresponding eigenstates |a,), |a,.), |@3), .... etc. of the energy are chosen as 
basic states, the Hamiltonian is diagonal, with elements aj, a3, a3, .... etc., and 
can be written as the sum of two diagonal matrices, one (H,) with elements 
G1, Go’, Ay’, 4x, ...., and the other (H,) with elements 0, @, — ay’, ag —a,', 0, pas 
where a, —a,/ and a;—a,' are small. Regarding H, as a small perturbation on /,,. 
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itis seen that the basic states appropriate to the matrix Hj are also| qa, ), | 4 ), layin 
Thus, the quantum mechanical description of a dipole with two nearly degenerate 
levels is equivalent to regarding the dipole as being slightly perturbed from a 
description in which the levels actually coincide. 

The reasoning which led to formulae (1), (2) and (3), was not dependent on 
the form of W; the assumption that W was of the form &,, -,W,,, was only made to 
find the moments due to inter-dipole couplings. In dealing with nearly degenerate 
levels, it is more convenient to regard the dipoles as having levels which actually 
-coincide, and to take the perturbation as a sum of internal couplings, suchas Hy, 
-and external couplings, such as dipolar couplings, and write W= Wine + Wort: 
The formulae (1), (2) and (3) are then directly applicable. The area, which does 
not depend on W, is given, as usual, by a sum of terms like that in (4), with one for 
-each transition. The first moment may be written as 


| exp(—£,/RT) Spur (PsW,.P)MP,M — P,Wiy,P,MPpM) 
Wo BT 4 exp(—E,/RT) Spur (P;Wex,P3MP,M — P,,WoxeP,MP 3M) |’ 


ext 


_and gives the first moment about v=(a,’—a,')/h. The first summation can be 
-evaluated as if there are no inter-dipole couplings, and therefore gives N times the 
moments of the transitions a,<—>a, and a,<—>a, of a single dipole about 
(ay —a,)/h. The second term is the first moment of a system whose dipoles have 
a degenerate energy level, and, as usual, is zero. It has not yet been necessary 
to give a value to a,’, except to require that a, — a,’ and a3 — a,’ are small, but if it is 
now chosen so that the first summation also vanishes, the total displacement 
vanishes, and the frequency v=(a,’ —a,)/h gives the centre of the absorption. 

In the same way, the second moment can be written as the sum of four terms, 
which may be characterized by the notation (W,4 Wint), (WextWext)s (WextWint) and 
(Win, Wext) where 

(WW) =u, xP (—E,/RT) Spur (PW, P,W_P,MP,M 


po 


—2P,W,P,MP,;W,P;M+P,W,P;W,P;MP,M), 
al’ » [MAGS 33 Bir ep BG 1B 2 


where W,, and W, can be W,,, and We. (WintWing) is the expression for the 
second moment which would be obtained if there were no couplings between the 
-dipoles, and is thus N times the second moment about (a,’ — a,)/h for the transitions 
-Qy~—> @, and a;<-> a, ofasingle dipole. (W..4W.,,) is the second moment which 
would be obtained if the dipoles had a degenerate energy level a,’, and is a case 
which has already been treated. When (W..,Wi4) and (Wiy,W..4) are evaluated 
itis found that they vanish. Taking the ratio of the second moment to the area, 
it is seen that the mean square width is equal to AW,? + AW,?, where AW,? is the 
mean square width of the system when there are no inter-dipole couplings, and 
AW,’ is the mean square width for a related system with a degenerate energy level 
ay, as discussed above. 


SG TRIALS IRUE ONE, SPIRO C IMU. 


In applying these results to the widths of the resonance lines given by para- 
magnetic ions, it is usually assumed that an ion can be completely specified by 
its electronic states, whereas, in fact, if there are couplings between the electronic 
system and the nucleus, a complete description will be a sum of products of nuclear 
‘and electronic states. If the nuclear—electronic coupling is zero, the ionic states 
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can be taken as |e) |n) where |e) is an electronic state and |n) is a nuclear state, so 
that each level, instead of being non-degenerate, as has usually been assumed, has 
actually a degeneracy arising from the nuclear states. If formulae (1) and (2) 
and (3) are evaluated, taking into account the nuclear states, and assuming that 
W is purely electronic in nature, it is readily seen that the only effect is to increase 
the degeneracies of P, and P, by factors which disappear when the ratios giving 
the mean displacements and mean square widths are taken. Thus the omission 
of the nuclear states does not affect the results of the previous work, provided that 
there is no coupling between the nucleus and the electronic motion. The assump- 
tion that this is so, is not, however, always justified. For example, the magnetic 
field at the nucleus of the Cu** ion, arising from the electronic system, is of the 
order of hundreds of thousands of gauss, and the interaction between its nuclear 
magnetic moment and the electronic system produces splittings of the nuclear 
levels of the order of 0:03 cm~. If the microwave resonances of the copper ions 
in a diluted crystal of copper potassium sulphate are observed, a hyperfine 
structure, with an overall separation of this order of magnitude, is found. It is 
necessary to use a diluted crystal since, in an undiluted crystal, the broadening 
due to inter-ionic couplings is also of this order, with the result that the lines in the 
hyperfine structure are broadened into each other, and are no longer resolved. 
Evidently then, in calculating the mean square width of the line in an undiluted 
crystal, it is incorrect to neglect the nuclear—electronic couplings. 

The necessary extensions of the theory in order to include the effects arising 
from hyperfine structures, are readily made. As for the case of nearly degenerate 
dipole levels, W is taken as a sum of external and internal couplings, W,,, (assumed 
to be entirely electronic), representing the couplings between the different ions, 
and W,,, representing the electronic—nuclear couplings which have given rise to the 
hyperfine structures of the levels. It is then found that the centre of the absorption 
line is at the same frequency as the centre of the hyperfine structure for the reson- 
ance of a single ion, and that the mean square width is equal to the sum of the mean 
square widths of the system with no electronic—nuclear couplings and the mean 
square width with no inter-ionic couplings. ‘Thus the results are completely 
analogous to those found for the case of overlapping absorptions caused by near 
degeneracies in the electronic levels. 


§10. DEPENDENCE ON THE SHAPE OF THE CRYSTAL AND ON 
TEMPERATURE 
The double summations of the form &;,|«¢|W,,|)|? which appear in the 
expressions for the main, temperature independent, part of the mean square width 
can be replaced by summations over a single index provided that the crystal is 
reasonably large. For, if a dipole z in the interior of the crystal is considered and 
the summation over 7 is performed, the series is so rapidly convergent that its 
sum depends, effectively, on only those dipoles in the immediate neighbourhood 
of iz. It will therefore be the same for all dipoles of the same type in the crystal, 
except for those within a few atomic distances of the surface, which, in general, 
form only a negligible proportion of the whole. ‘The double sum is thus N times 
the single sum. It will be noted that for dipolar interactions, |{|W,;| >|? falls 
off as the inverse sixth power of the distance between the dipoles; for exchange it 
falls off exponentially. ‘The convergence is therefore very rapid, and the final 
sum does not depend on the shape of the crystal. 
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The terms of order 1/RT which appear in both the mean displacement and. 


the mean square width of the line, contain summations of type Xi 5 «| Wil D- 
For exchange interactions such sums are rapidly convergent, but for dipolar 


interactions the convergence is no longer so rapid, and it is no longer possible 


to treat the double sum as N times a single sum. The contributions of the 
surface dipoles are, in fact, appreciable in the interior of the crystal, and the total 
sum will depend on the shape of the boundary. 


In deriving the general formulae for the moments it was assumed that the 
temperature was so high that it was a good approximation to assume that. 
exp {—a|W]|a)/RT} could be replaced by unity. In discussing terms in 1/R7” 


this approximation is no longer valid. A more accurate procedure would have 
been to expand exp { — (a| W|a)/RT} in powers of 1/RT, retaining all the terms. 


When this is done itis still possible to write the moments in invariant forms though 


they area good deal more complicated. ‘The expression for the area is quoted as an. 
example. 


Area 


Spur P,WP,MP,M 


x, pexp (— E,/RT) { Spur P,MP,M — = 


1 
ots ORT? Spur P,WP,WP,MP,M.. 5 +. 
In principle any sum of the form 


yexp(—E,/RT) Spur Z, 


where Z contains P,, P;, W and M in any combination, can be reduced to a 


summation of matrix elements of the couplings between pairs of dipoles. In 
practice only the simplest expressions can be evaluated without prohibitive labour. 
A general type of sum will have the form 


X15... Weld Wad. D- 


If there is a sufhx occurring once only, summation over that suffix will lead to a 
shape-dependent sum. If the suffix is repeated, summation over it will give a 
sum whichisindependent of shape. The approximation exp {— (a|W|a)/RT}—1 
is therefore likely to be a good approximation provided that a summation like 
=,¢| W,;|) remains small compared with RT, for this ensures that each term in 
the series for the moments is smaller than the preceding term, or that the series are 
rapidly convergent. ‘The first term is then a good approximation to the actual 
value of the series. 

On a classical model, a summation like &;W,; measures the energy of the th 
dipole in the field of all the other dipoles. This internal field varies as 1/RT and 
is small at temperatures well above the Curie temperature. It is also shape- 
dependent and is usually described, for a uniformly magnetized body, in terms of a 
demagnetizing factor. ‘Thus it is reasonable to assume that the formulae obtained 
for the moments are good approximations provided that the temperature is well 
above the Curie temperature. Near the Curie temperature the line shape and its 
displacement are likely to depend on the temperature and on the shape of the 
crystal. 

In a first approximation the effect may be regarded macroscopically as arising 
from the magnetic field of the moment induced by the external field (which is 
proportional to 1/RT). The steady component (i.e. the demagnetizing field), 
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superimposed on the external field, changes the effective steady field on a dipole. 
The demagnetizing field corresponding to the oscillating component depends on 
frequency in a resonant way, and taken with the oscillating applied field causes a 
shift in the resonance frequency. 


§11. FREE SPINS 

The simplest application of these formulae is to the case of dipole-dipole 
coupling between free spins, such as occurs in nuclear resonance experiments. 
It will be assumed that the system consists of a mixture of spins, some with spin S, 
and g-factors g, and others with spin S, and g-factors g,. In the absence of any 
couplings, the microwave spectrum will consist of two lines at frequencies g,BH/h 
and g,8H/h where f is the Bohr magneton and H is the external magnetic field. 
Assuming that all the spins of the same type are equivalent, and using the forms of 
equations (7) and (8) appropriate to the cases in which several transitions of the 
same spin occur at the same frequency, the mean square width of the line at 


o,BH/h is 


3 gi*BAS(S, He) ee COS Oy \" a 1 gy'go?B*So(S2 + Lye ef Ecos Cie! 
4 h? wd : 13 h? fe ate 


Here, 7; is the distance between spins | andj, 6,, is the angle between the radius 
vector 7,; and the direction of the external p.c. magnetic field, spin 1 is a typical 
spin of the first kind, 7 ranges over its similar and k over its dissimilar surrounding 
spins. This formula is due to Bloembergen et al. (1948). It then appears that, 
apart from the factors g,2,S,(S, +1) and g,2,S,(.S, +1), similar ions give 9/4 times 
more contribution to the square width than do dissimilar spins, a difference due to 
the resonance effect. 


at 


12 EXCHANGE INTERACTIONS 
An interaction of particular interest is that which arises from electron exchange 
betweenions. For the present discussion we confine ourselves to ions whose spin 
is 4. The exchange energy of two ions, z and 7, can then be written as: 


Wi; = 2J 5; S, . S; 


(where J,, is independent of the spins). 

If all the ions are identical and this form of interaction is inserted into the 
appropriate form of (5) the second moment is found to be zero. ‘This implies that 
exchange interactions alone give rise to no broadening effects. ‘That this is so can 
also be seen from the following general argument. The exchange energy of the 
system 

We = Lj. 
commutes with all symmetrical functions of the ionic spin variables, and hence 
with the unperturbed Hamiltonian, #%,) (the Hamiltonian of an uncoupled 
system). The total Hamiltonian is therefore the sum of two commuting terms and 
its eigenvalues are the sums of the eigenvalues of # and W,,. The correspond- 
ing eigenstates are simultaneously eigenstates of #) and W.,. The probability 
of a transition being induced by a weak transverse oscillating magnetic field is 
determined by the matrix elements of the transverse magnetic moment. ‘This is 
also a symmetric function of the spin variables and therefore commutes with W,.. 
Thus it has no non-vanishing matrix elements connecting states of different 
exchange energies, or, in other words, the selection rules are AW,,.=0, 
AH, = +gBH, giving rise to one sharp line. 
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If the ions are not all identical the above discussion is no longer applicable. 
There are two cases of practical interest : (i) the resonance frequencies of the ions 
are not all the same, (ii) the ions have the same resonance frequency but are 
magnetically anisotropic with their principal axes not allthe same. We consider 
these two cases in some deta‘l using a system in which there are two types of ions. 
The effect of exchanye in the first case is to give a broadening proportional to the 
‘exchange frequency’, J/h. Just as an isotropic ion can be regarded as precessing 
about the direction of the magnetic field, so an anisotropic spin precesses about an 
axis which, in general, is different in direction from that of the magnetic field. 
In the second case the two types of ions have the same precessional frequencies but 
are precessing about different axes. The broadening is then proportional to J/h 
and to sin? 46, where 6 is the angle between the two precessional axes. 

In the first case we suppose that there are two types of ions, one type precessing 
with frequency v about the z-axis and the other with frequency v’ about an axis 
inclined at an angle @ to Oz. The relevant matrix elements of W;; are 


(1, 1] Wil 1,1) = (2, 2] Wy] 2,2) 
== 2,1 | W,;|2,1)= — 1,2] W,,|1,2) =4,,; cos 0, 


where ionszand/are of different types. Inserting these values into the appropriate 
form of (5), the mean square width of the line at v is found to be 


Av? = cos?6 D,J?,/h?, 


where ion 1 is a typical ion of the first type and the summation is to be taken over 
all ions of the second type. This analysis is of course only valid if the two lines are 
distinct, i.e. provided that J cos @ is sufficiently smaller than A|v—v’|. When J 
is larger the two lines fuse and the width can be calculated by another method 
(Pryce 1948). 

In the second case the frequencies of the two types of ions are the same and the 
appropriate form of (5) now contains the matrix element 


<1,2| W,,|2, 1) =J;,; cos? 46. 
In this case 


(Av?) =sint $0. UJ, 2/h?. 

Reverting now to the case of identical ions, it has been shown by van Vleck 
(1948 a) and Gorter (1947) that the combined effect of dipolar and exchange inter- 
actions is to make the resonance narrower than it would have been with dipolar 
coupling alone. ‘The narrowing effect of exchange interactions can now be 
compared with the narrowing effect in nuclear magnetic resonance in liquids 
arising from the rapid molecular motions. Their argument is that the second 
moment is unchanged by the presence of exchange but that the fourth moment is 
increased. ‘I'his can only happen if the absorption consists of a peak which is 
sharpened and wings which are simultaneously extended by the exchange. 

‘The argument can be extended to indicate that when the exchange energy 
is much larger than the dipolar interaction energy, the absorption is nearly all 
concentrated in a narrow peak, a small fraction being spread out into the wings 
over a frequency range of order J/h. 

To see this, an interaction energy 


W=dW,, +pWaip, 


is considered, where A and p are parameters which will eventually be put equal to 


The Theory of Magnetic Resonance—-Line Widths in Crystals 51 


unity. It is readily seen that the 2nth moment will be a polynomial in A and pu of 
the form 


h2"(Av)2” = Ad2n-2),2 + BYn-3y3 4... + Cu2”, 
where A, B, C etc. depend on z and on the coefficients in the exchange and dipolar 
energies. The termin A?" is absent, for if  =0 it is known that (Av)?” = 0 (exchange 
only). ‘The term in A?”“1, is absent because the polynomial must be non-negative 
for all values of A and y.. If the exchange effects predominate over the magnetic, 
the polynomial is dominated by the term in \2”-22, whose coefficient will contain 
a factor J?”-*62H?, where/J, 8, and H, are typical values of the exchange interaction, 


Magnetic moment, and internal magnetic field. Setting \ and w equal to unity 
gives 


h2n(Ay)2" = a, J?" BH, 
where a,, is some numerical factor of order of magnitude near unity. This is 
consistent with a resonance in which nearly all the intensity is in the central peak, 
which does not contribute appreciably to moments above the second, together with 
an intensity of order 8?H;?/.J? spread over a range of frequency of order J/h. 
This argument throws no light on the width of the central peak. 


§13. ORIENTATION OF THE OSCILLATING MAGNETIC FIELD 

When there are no coincident frequencies it is readily seen that the transition 
probabilities |{1]|2)|? drop out of the mean square width; the mean square 
width does not then depend on the orientation of the oscillating magnetic field. 
If, however, there are coincident transition frequencies, they do not drop out, and 
further they occur as actual matrix elements rather than as square moduli of matrix 
elements. Although the oscillating magnetic field may be at right angles to the 
steady magnetic field, so that the transition probabilities are the same for all 
positions of the oscillating field, the actual matrix elements will change as the 
direction of the oscillating field is changed in the plane perpendicular to the 
steady field. This means that the mean square width will depend on the orienta- 
tion of the oscillating field. 

§14. CONCLUSION 

The purpose of this paper has been to give the principles and general results of 
the method rather than apply it to specific cases. There are a number of cases for 
which the widths have been measured experimentally and the comparison of the 
theoretical predictions with the experimental widths will be given later. 

The subject of this paper has recently also been treated in an important 
paper by Van Vleck (1948b), which covers much the same ground. Our 
method of approximation differs in detail from Van Vleck’s, and in some direc- 
tions our results are more general than his, while conversely, in other directions, 
notably the explicit calculation of the fourth moment, his results are more 
complete than ours. 
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ABSTRACT. A theory is developed for the capture of particles into synchrotron orbits 
when the radio-frequency accelerating voltage rises from zero to its maximum value in a 
finite time. It is shown that if this time of rise is sufficiently long, all particles originally 
occupying a band of width 1/./2 times the maximum width of the fina] stable region for 
synchrotron phase oscillations wil! be captured, irrespective of the initial phase at which they 
enter the accelerating gap. ‘This conclusion would seem to be in agreement with the 
observed characteristics of existing electron synchrotrons. 


Sail JUNAERKONDIOC TION 
ARTICLE acceleration employing radio-frequency accelerating voltages 
combined with a magnetic field to define circular orbits was first achieved 
in the form of the cyclotron. In recent years the principle has been extended 
to the acceleration of electrons (electron synchrotron) and of heavy particles to 
higher energies (synchrocyclotron and proton synchrotron). 

Because of relativistic effects it is in general necessary to vary two of the 
parameters H, r and w (magnetic field, orbit radius, and angular radio frequency) 
during the course of the acceleration. In the synchrocyclotron, r and w are 
varied, except in the case of non-relativistic particles (cyclotron), when w may be 
kept constant, while in the synchrotron, H and w are varied, the exception this 
time being for extreme relativistic particles (electron synchrotron), when only 
H need vary. The particular advantage of the synchrotron is that with the con- 
stant equilibrium orbit radius the magnetic field need only be established in a 
restricted region about this orbit, resulting in a considerable saving for high 
energy accelerators. 

Inthe electron synchrotrons operatiny at the present time it has been convenient 
to inject the particles at velocities considerably less than the velocity of light, 
and to accelerate them to extreme relativistic energies (of the order of 2 Mev.) by 
a betatron mechanism. As this betatron action slowly diminishes (by magnetic 
saturation) the betatron orbit contracts, and when its radius matches the synchro- 
tron equilibrium radius, the R.F. accelerating voltage is applied to a resonator 
surrounding the orbit. As the betatron action further diminishes, the deficit 
is made up by the R.F. voltage. 

Not all of the particles initially present can in general be captured into stable 
synchrotron orbits, and this paper discusses the effect on the capture efficiency of 
the time of rise of the accelerating voltage, a lower limit to which is in practice 
set by the selectivity of the resonator. 

In general it is found that if this time of rise is above a certain limit, considerably 
more particles may be captured than would otherwise be the case, which in fact 
would account for the capture efficiencies approaching 100° observed in 
all existing electron synchrotrons. This is of greatest importance for schemes 
involving repeated capture, i.e. handing over of the particles from one resonator 


or R.F. harmonic to another, as is necessary, for instance, in the accelerator proposed 
by Kaiser and Tuck (1948). 
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In the case of the proton synchrotron, the critical point is likely to be in the 
initial stages of acceptance and acceleration. Some relief from these critical 
conditions seems to be possible by selecting the time of injection and the rate of 
rise of the R.F. voltage. 


2. THEORY OF PARTICLE CAPTURE 

The theory of particle stability and bunching in the synchrotron has been 
developed by several authors (Bohm and Foldy 1946, Frank 1946, Dennison and 
Berlin 1946). ‘The most useful treatment, giving a clear physical picture of the 
process, is probably that of Bohm and Foldy, who consider a particle of charge e 
moving with a velocity § times the velocity of light in the r, 6 plane of a cylindrical 
polar coordinate system. In this plane the magnetic field has only a x component 
defined as a function of r by 

d(In #) 
es apie eee (1) 

The phase at which the particle enters the accelerating gap (relative to the R.F. 
phase) is 


pails dh MAR | tee (2) 


while a particle which always travels on the equilibrium orbit must enter at the 
equilibrium phase given by 
Sip OL Ny elt ee eee (3) 
where , is the angular velocity of the equilibrium particle (and the angular 
frequency of the r.F.), V is the peak r.F. accelerating voltage, and wv is that part 
of the total accelerating voltage per turn required to maintain a particle on the 
equilibrium orbit not supplied by betatron induction. 
For a particle to be successfully accelerated its phase must execute damped 

oscillations about ¢,, and the equation defining these oscillations is 

d { 2rk, 

F (sae) Hey Sit O= VSI Oa Sad eens (4) 
where E£ is the total energy of the particle (including the rest energy) and 


n 
eee fella): 
the subscript, s, everywhere referring to the equilibrium values of the quantities. 
Equation (4) is analogous to that defining the equation of motion of a pendulum 
of moment of inertia J =27E,/ew,? K, restoring torque V sing, and constant bias 
torque v=Vsin¢,. Since £, will in general remain substantially constant over 
a period of the phase oscillation, we can, for the quasi-static case, rewrite (4) as 


$b +kV sing =kV sind, Sy tes (5) 
with k= ew,?K/27E,, which, on integration, gives 
$?/2kV =cosh—cosdm+(b—bm) SING ve eee (6) 


where ¢,, is the maximum phase reached during the oscillation. 
Associated with the phase oscillation is a radial oscillation defined by (6) and 
1 ee a ee oe ee (7) 
r, PB (1—n)Ko, 
where « =7-—7,. 
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We can rewrite (5) as 


ees 


_ oU(¢) 


Sous 8) 
kV of ( 
where U(¢) is the potential function defined by 
U(¢)= —(cosd+¢sing,), = snes (9) 


and is plotted for various values of 4, in Figure 1. The maxima in these curves 
occur at 6=7-—4¢, and the minima at $=4,. 
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Figure 1. The potential function, U(¢), for various values of ds. 


The limiting condition for stable oscillation is found by putting ¢,,=7—¢, in 


(6), and we find that stable orbits only occur for é and d within the region of phase 
space enclosed by 


b= + (2KV[(cosd +cos $,)—(7——4,) sin gJ}4, (10) 
which equation actually defines the shape of the bunch of particles, taking ¢ as 
azimuth relative to the bunch centre. 

The limits of the stable region defined by (10) are plotted in Figure 2 for various 
values of $,, while Figure 3 gives an idea of the motion of the particles in phase 
space for both zero and non-zero values of 4. 

The maximum half-width of the stable region is found by putting 4=4, in 
(10), giving 

ou =(4kVO)?, re Ae, 
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where ® = cos ¢;—[(7/2)—¢4,] sin ¢,, and ¢y is the maximum possible phase 
velocity within the stable region. Figure 4 gives ©! as a function of V/v. 


$(4RV)-4 


10 


Figure 2. Envelopes of the region within which 
stable phase oscillations are possible for 
various values of ¢s. (b) 


Figure 3. Particle trajectories in phase space 
for (a) ds=0°, (6) ds=30°. The dotted 
curves give the boundaries of the stable 
regions. 


Figure 4. ©? =dy,/(4kV)* as a function of V/v. 


As well as the radial oscillations associated with the phase oscillations there 
are also rapid free oscillations (betatron oscillations) which, however, have been 
neglected in arriving at equation (4), since the coupling between the two motions 
is in general negligible. Because of this, in the whole of the following the free 
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oscillations will be neglected and the particles will be considered to be travelling 
on their instantaneous orbits, even when pure betaton acceleration 1s considered ; 
thus, by the width of the betatron beam will be meant the spread of betatron 
instantaneous orbits. 


(i) Instantaneous Capture 


The above considerations enable us to determine the transition from betatron 
to synchrotron acceleration when we instantaneously* establish the accelerating 
voltage. This case has been worked out by Goward (1947). Prior to transition 
we have a betatron beam of width ¢, and subsequently a stable region of half- 
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Figure 5, Illustrates particle capture with an instantaneously established r.F. voltage, for 
the ideal case ¢ds=0°. 


(a) Betatron beam of half-width ¢. 
(b) y= <€o, 64% capture efficiency. 
(Cc) €y;=5€o, 94% capture efficiency. 


width ey, and only those particles will be trapped which at the instant of transition 
lie within this stable region. At this instant the betatron acceleration will generally 
be almost adequate to maintain the acceleration, so that 4, can be taken as zero 
and the process is as depicted in Figure 5. It is clear that for efficient capture we 
need ey >< 9, while for ey, ~€) we only get about 60° capture. Subsequent to 
transition the particles will oscillate about the equilibrium radius with radial 


* Actually equation (4) is derived by making a Fourier decomposition of the accelerating gap 
voltage into a series of travelling waves and selecting the component in step with the particle motion, 
which is the only effective term over a few revolutions. This, however, represents an averaging 


process, so that the term ‘instantaneous’ implies ‘in the shortest time for which (4) will be valid >— 
in fact, in the time of a few revolutions. 
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amplitudes from zero to ey, and as the betatron flux saturates further, f, increases 
and the stable region contracts, both in azimuth and in radial width, causing a 
fraction of the particles to be lost. For an increase in ¢, of from zero to 30° the 
fraction lost may amount to 40%. 


(11) Finite Time of Rise of the R.F. Voltage 


An examination of the parameters of a number of existing electron synchrotrons 
shows that in no case does the R.F. voltage rise instantaneously in the sense 
indicated above, and that in a conventional machine it may be difficult to satisfy 
this condition. 

In the following, an attempt is made to solve the phase equation (5) for the case 
in which V rises from zero to a value V, in a time t,. While a general solution 
would seem impossible, we can obtain an approximate solution, subject to 
restrictions to be placed on 4, as follows: 

From (5) we have 

2 
a TE 2RIasine = 2ko, Pe Oh ee (12) 


and integrating, 


2 = 4,2 +2kV cos db + 2ku( — bo) — iE eh 


where at time t=0, $=4,, ¢ $y, V=0. 
The problem now is to determine the range of values of ¢) and ¢y which a 
particle may have in order that it may be captured by the time ¢=4,. Capture 


actually occurs when d becomes bounded, i.e. the motion becomes oscillatory. 


Inspection of (13) shows that it consists of a term bo? + 2hv(b — b); which is a 
linear function of ¢ plus an oscillating function of ¢ consisting of two terms, the 
predominant one being 2kV cos¢ if V is sufficiently small, with amplitude 
increasing with time. ‘Thus $2 has a series of minima, and capture may be said 
to occur at that minimum when 2 just becomes zero, since after this it will be 


bounded and ¢ will oscillate about ¢,. 
From (12) these minima occur at values of ¢ given by 


sing =sin ¢, | rest (14) 


cos ¢ = cos ¢, 
or 

Ge (20 li eyo 0 ee eee (15) 
where 

Cp AUS il ee pet Pe cee 

Substitution of (15) (with the appropriate value of q), plus the simultaneous 
condition ¢?=0, into (13) then gives the final result that a particle with initial 
coordinates ¢,) and ¢y will be trapped when the voltage V satisfies 


0 =d,2—2kV cos d, + 2kv[(2q+1)7—-b,-bo]-R, +e (16) 
where 


eRVpOS OY 
$ db : 


ae 2kVc 
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It may seem that for particles which are not trapped the above conditions ae 
satisfied at the point where they cross the equilibrium orbit and ¢ changes sign. 
(see Figure 3). In such cases, however, d?/d¢ is non-zero and hence (14) is not 
satisfied. 

Equation (15) is not satisfactory since we do not know the appropriate value 
of q (i.e. the number of phase revolutions prior to capture); however, one thing 


that is clear is that by making V sufficiently small (i.e. ¢, sufficiently long), R can 
be made negligible compared with the term 2kV cos¢,._ In what follows, then, 
R will be treated as negligible and later an estimate will be made of its magnitude 
for various initial conditions, and the restrictions to be placed on ¢, determined. 

Returning to (13) and writing ¢, =(2v + 3)z (v integral, positive or negative), 
we get, neglecting the integral term, 


$, = go +2ku(d, — bo); 
and thus $2 =4,2+2kV cosp+2ku(p—4,), eee (17) 
which with (15) and de? =() gives the trapping condition 


0 =¢4,2—2kV cos 4, + 2kv (5 -$,). 


For final conditions V=V,, ¢,=¢,, all particles will be trapped which le 
within the limits of by defined by 


$,2 =2kV, cos f,, — 2kv G ~$) 2h Os (18) 


Now be is obtained by putting ¢=¢, in 
d= $5? + 2ho($ — do), 


corresponding to $ =k, ites $ = by + hut, and since the time corresponding to 
d=, is approximately ¢, (for V sufficiently small) we can write 


, S $6 = kvt,. 
Thus all those particles will be trapped which have bo within the range 
y= — vt, + (2kV,®,)}. 


The important point of this result is that it does not contain ¢). Further, 
by comparison with (11), the range of radii corresponding to this range of dy is. 
just 1/./2 times the width of the final stable region, and all particles within this. 
range will be trapped, irrespective of the initial phases at which they enter the 


resonator gap. This means that for 100°, capture the resonator voltage only 
needs to be sufficient to make 


Ey SA 26). ae) Ate eee (20) 
instead of >, for instantaneous capture. 

The term —kvt, in (19) represents the average orbit contraction during the 
time ¢, of establishing the r.F. voltage. This will set an upper limit (which will 
be trivial for most electron synchrotrons) to the time of rise, since this amount 
of contraction must be possible within the orbit space. 
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We can now determine the equation of motion, immediately subsequent to 
capture of a particle whose initial phase velocity satisfies (19). At the instant of 
capture we have 


$2 =2kV, E due é = $a) sin $a | 


$2=4,2 + 2kV cos $ +2hu(b —F,). 
Substituting for b. and putting ¢’ =¢ —2gz7, we get 


¢? =2kV, [cos $,, +cos d’ —(7-b —¢') Sindy], ~ ....-. (21) 
which reveals by comparison with (10) that such particles, as we would expect, 
oscillate with the maximum stable amplitude, i.e. on the perimeter of the stable 
region defined by (10). 

An estimation of the conditions under which the integral term R (equation (16)) 
may be neglected, and thus (19) is valid, is given in the Appendix. The most 
important result is that if V,/vS1, R is negligible when 


CGR) ae leo ee es (22) 

and further, that if DARD MCR 0 )8 rae = ee PP Fr ea (23) 
both R and kvt, are negligible and (19) becomes 

(Nos C155) vat Bian a) Rte ae (24) 


Equations (23) and (24) are equivalent to the statement that if the time of rise of 
the resonator voltage is long compared with 1/47 times the period of phase 
revolution of the extreme particles prior to capture, and yet sufficiently short that 
the orbit contraction during this time may be neglected, then all particles with 
initial phase velocity within the limits of (24) will be captured, irrespective of 
their initial phases. 

Equation (22) will be the condition for most electron accelerators, since at the 
time of transition v in general will be small and V,/v>1. If this is not so 
(as may be the case for accelerators in which initial betatron acceleration is not. 
employed) the appropriate condition may be found from the treatment in the 
Appendix. 

(iii) Suppression of the Resonator Voltage 


In a similar fashion to the above, we can show that if the R.F. voltage falls. 
from V, to zero in a time ¢, satisfying (22) (or the appropriate criterion given in 
the Appendix if V,/v+1) the particles will be spread into a beam, uniformly 
distributed in phase, of width 1/1/2 times the width of the stable region prior 
to suppression. 

On the other hand, if the R.F. were instantaneously suppressed, the particle 
bunch would retain its shape immediately subsequent to the suppression, but would 
then rapidly become uniform in phase, its width being equal to that of the stable 
region prior to suppression. 


(iv) Saturation of the Betatron Flux and Damping 


Subsequent to capture of the particles into stable orbits there are two major 
effects. The first is the contraction of the stable region due to the magnetic 
saturation and the second is the phase oscillation damping as the magnetic field 
(and particle energy) increases. 
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The phase oscillation damping for the case 6,0 is given by 
VR 
B¥(byn/2) sin ($y/2) & keg = © igus (25) 


where B is the complete elliptic integral defined by 


te cos? u du 
BG n2)= |. A asint Gamba 
Equation (25), together with 
bm/out = €m/€m =sin (¢m/2), Maa Seas (26) 


defines the damping of the associated radial oscillations relative to the stable 
region (where e,, is the radial amplitude and €y, is the maximum half-width of the 


stable region). 
For small oscillations B ~7/4=constant, and we get 


and 


woe2K |* 
Em/€u we | a | a eS ie ratrioa. (28) 


The damping corresponding to increasing particle energy is plotted in Figure 6, 
in which the energy has been normalized to unity at ¢,,/ey;=1. The broken line 
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Figure 7. Illustrates the effect of betatron saturation 
assuming v/V to increase linearly from 0 to 2, 


Figure 6. Phase oscillation damping with while Es increases from 1 to 2 (relative units). The 
increasing particle energy, Es. The curve (a) gives &?=¢,,/ (4kV)?, which corresponds 
conde corresponds to damping to the limit of the stable region, while the dotted 

Ss 4. 


curve (b) gives bml(4kV)*, corresponding to the 
amplitude of oscillation of the particle which is 
just retained during the process. 


‘corresponds to «,/ey <E,*, and inspection shows that the damping for large 
amplitudes is less than this, the difference for many purposes not being appreciable. 

The saturation of the betatron flux in general takes place slowly, thus v (and 
-p,) will increase from zero to some maximum value causing a contraction of the 
‘stable region (see Figure 2). The purpose of this section is to determine the 
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area of the initial stable region from which particles may be lost, and it is shown 
that if V, is made somewhat larger than needed, just to trap all of the particles, 
this area will contain no particles, and thus there will be no loss (ignoring possible 
space charge effects). 

Since ¢, will generally be negligible at the time of transition, we really wish 
to know the effect of slowly increasing ¢, from zero to a maximum value Pap? 
Consider for instance the case for which 4, has initial and final values 0° and 30° 
respectively, which is probably about as large a change as would be experienced 
in a conventional electron synchrotron, An examination of the potential 
curves of Figure 1 shows that the particle which is just lost during such a change 
is one which initially corresponded to a value of U(¢,,) within the shaded area. 
It is clear that we will not be seriously in error in taking the upper limit as 
representing the initial level for this particle, i.e. in assuming that 


U(¢n) — U(bs) =bn2/2RV 


remains unchanged during the change of shape of the potential well (which is in 
fact the case for small oscillations). 

Taking into account the phase damping during the period of magnetic satura- 
tion, which we will assume to be of the order of £* (even though large angle 
oscillations are involved), we see that [¢,,/(2kV)]E remains constant during the 
change in ¢.. If ¢, reaches a maximum value ¢,,, all particles with peak phase 
velocities greater than 


bm= (2) CL a es (29) 


will be lost, where the subscript p refers to the value when 4,=4,,, all other 
quantities being evaluated at transition. 
Translating phase velocity into radius, (29) becomes 


1 E, D t 
Em — ®,, (=) ew stsitotls) (ei = G 0) 


(remembering, of course, that ¢,,>€,,). 
The process is illustrated in Figure 7, which depicts the case where v/V 
increases linearly from 0 to 2, while in the same period the particle energy is 


doubled. sy 
Let us now consider the case of a betatron beam of half-width «) which is 


captured by a slowly rising resonator voltage V. ‘The extreme particles of phase 
velocity do will be captured when V = Vo, where ¢)=(2RV))?. As V increases to 
its final value V,, the phase oscillations will be damped approximately as V-*, and 
at the same time on will increase as V?, so that bm (and «,,) will increase as V+. 
Thus in order that the extreme particles are not lost during the period of saturation 
we must have 
De Vga hg Ba | tee (31) 

The relation between the width of the betatron beam and the width of the stable 
region immediately subsequent to transition, in order that all of the electrons be 


retained, is now 
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(where €);=¢€ yy when V =V,), which replaces the condition (20) in cases in which 
there is a significant change in ¢, subsequent to transition. We must again 
remember that ¢) cannot be greater than ey,/+/2, ie. if (#,/E,)'®,>1, the 
appropriate relation, is (20). 

The above result can only be regarded as approximate because of the assumption 
that the damping is always the same as for small angle oscillations, but its signi- 
ficance is in showing that if the resonator voltage rises slowly to a value somewhat 
greater than is needed just to capture the betatron beam at the time of transition, 
there will be no subsequent loss of particles during the magnetic saturation. 

The condition that the theory of capture be valid is, of course, 


E(ZEV GS Se eee (33) 
where t=f, when V =V , which is equivalent to 


HOV >a) a (34) 


§3. APPLICATION TO ACTUAL ACCELERATORS 


An analysis of a number of existing electron synchrotrons has shown that in 
all cases the time of rise of resonator voltage satisfies (33), which is in general 
agreement with the fact that the capture efficiency in all cases approaches 100%. 


(i) 70 mev. Synchrotron 


An excellent example is the 70 Mev. machine described by Elder et al. (1947), 
who state that increasing V,; beyond 2v, does not appreciably increase the x-ray 
output. 

The parameters of interest are: v,=170 volts and thus V,=340 volts, 
w,/27 = 163 Mc/s., n=0-75, and at transition E,=2-5Mev., giving (2kV,)?= 
2x10%sec-!. The value of V/v at transition is not given; however, the half- 
width of the curve of output against time of switching on the R.F. is of the order of 
50 psec., and assuming that this is the order of the time of the orbit radius, in the 
absence of the accelerating voltage, to contract by an amount equal to the width 
of the betatron beam, we get V,/vu~ 200 at transition, and thus ¢, at this instant 
may be taken as zero. 

The ratio E/E, is not given, but as a rough estimate we shall put it equal 
to 2, and we should not be greatly in error in calculating «9 since this ratio only 
appears in equation (32) to the power one-half. 

Putting V,/v,=2 and £,,/E,=2 in (31) we get V,/V)=4. For the capture 
theory to be valid we need, using (34), ¢, > 0-2 wsec., which should be satisfied by 
t,~2psec. or longer. In fact, the above authors quote 100% capture for t, 
between 2 and 20 sec., while for an instantaneously established R.F. voltage we 
would expect an overall transition efficiency of the order of 50%. 

From (7), (11) and (32) we find that the half-width of the betatron beam prior 
to transition is _g~0-2cm. Since the spread of instantaneous orbits during the 
betatron acceleration decreases approximately as H~ (i.e. as BE) (Kerst and 
Serber 1941), the above figure corresponds to a spread of orbits at injection 
(injection voltage of 40 kv.) of the order of 4cm. It is interesting to compare this 
figure with the width of orbit tube, namely 9cm., indicating that the spread of 
orbits is of the order of one-half this width, a result which might be expected in 
the absence of effects of space charge, azimuthal field irregularities, etc. 
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In conventional electron synchrotons, irrespective of the final energy, there is 
not a great deal of variation in the betatron injection voltage, transition energy, and 
the ratio of orbit tube width to orbit radius, the main variables being the orbit 
radius and the peak voltage per turn needed to maintain the acceleration, Ue 

Taking the 70 Mev. machine to be typical, (33) may be reduced to 


J STALD y g e l (35) 


‘This means that the condition for the capture theory to be valid is that the R.F. 
voltage rises to Vy in a time which is long compared with ten times the R.F. period, 
a condition which will in general be satisfied in electron accelerators simply 
because of the high selectivity of the resonator, quite apart from limitations set 
by the associated electronic equipment. 


(ii) Heavy Particle Accelerators 


The capture theory is directly applicable to the heavy particle synchrotron, 
especially if initial betatron acceleration is contemplated in order to reduce the 
range of frequency modulation and to take advantage of the damping and con- 
traction of the instantaneous orbits during the betatron phase. In this case a 
suitable choice of the R.F. voltage and its time of rise will give 100° capture and 
can compensate for the phase antidamping with increasing w, in much the same 
fashion as the compensation for betatron flux saturation in the electron synchrotron. 

If injection is to take place directly, without the initial betatron phase, using 
an injector placed outside the equilibrium orbit, the situation is more complicated. 
The oscillation amplitude subsequent to injection is the sum of the free and 
synchrotron amplitudes, and since the damping of these free oscillations in one 
phase oscillation period will be small, an appreciable phase acceptance angle can 
only be achieved by first injecting the particles and then establishing the R.F. 
voltage, otherwise only particles injected near the equilibrium phase will be able 
to clear the gun structure. We could, for instance, inject until the orbit of the 
particle executing no free oscillations contracts to the equilibrium orbit, when the 
spread of instantaneous orbits will be a small fraction of the width of the orbit 
tube (due to the relatively small betatron induction), and the free oscillation 
amplitudes will extend from zero up to amplitudes just clearing the injector. 
If we then instantaneously establish the r.F. voltage (which, in contrast to electron 
accelerators, is easy to accomplish) the fraction of the particles accepted depends 
on the maximum phase acceptance angle and the radial width of the stable region. 
Depending on the machine parameters, there will be an optimum accelerating 
voltage corresponding to a stable region somewhat narrower than the available 
orbit space. For the Birmingham synchrotron the optimum acceptance fraction 
seems to be about one-third. 

On the other hand, with a slowly rising R.F. voltage we can accept particles of 
any initial phase and yet have a final stable region which is narrower than would 
otherwise be necessary, thus accepting particles with relatively greater free 
oscillation amplitudes than in the previous case. The drawback here is that since 
V,/v will be fairly small at the time of capture, the average orbit contraction during 
the time of rise may need to be a substantial fraction of the width of the orbit space. 
We may thus need to displace the equilibrium orbit towards the inner wall of the 
orbit tube and establish the radio frequency when the displacement of the 
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instantaneous orbits from the equilibrium orbit equals this desired orbit contrac- 
tion. The gain in injection efficiency may not be great, but in this case the particle 
oscillations may have a considerable clearance from the injector, allowing a 
relaxing of the tolerance on the radio frequency immediately subsequent to 
injection. 

The capture theory is also applicable if we desire to change R.F. sources. 
during the acceleration (e.g. from a varying to a fixed frequency source for 
extraction purposes). By slowly suppressing one source and switching on the 
other in the same fashion at a frequency corresponding to a suitable reduced | 
equilibrium radius, the difficult phasing problem disappears, leaving the relatively | 
simple matter of starting off the two events at the right moment. 


APPENDIX 
Circumstances under which the Integral R (Equation (16)) may be neglected 


First let us consider the simple case when v is sufficiently small that V,/v>1 
and at the same time the average orbit contraction duriny the time f¢, is negligible, 
Le, froma(l Oe 


kot, <(2kV)! 
or L(ORV)§<2V 0. ee (36) 
Equations (17) and (19) now become 
Gi= 2 P2kV cosh ene eeeren (37) 
do = £(2kV 8. ee ee (38) 


the fractional error in (38) being approximately R/4RkV,. 
Substituting (37) in the expression for R we get 


‘ Tt qd Te 
R=2nv | | Seah 2 =o JsC0S ae if 
J. do v=1/—2(¢,7+ 2RV, cos 6)3 
where 0=¢—¢,+4+7/2. 


Assuming V to be constant at V, over the periods 27 of 0, all of the elliptic 
integrals comprising the series drop out and we are left with 


2kV . 
R-=- —— SiN Pp. elfe) eral (39) 
$o 
Putting dy =(2RV,), the condition for (38) to be valid becomes 

iL(ZRV,)>3.. > | Re (40) 

Combining (36) and (40) we have 
2Vi0>t(2RV4)) 5, ee (41) 

which is equation (23) of the text. 


From (39) we can actually make a second approximation to the limits of bo 
within which particles are trapped, which become, instead of (38) 


b 


ee 74 | 7 = Sino 
bo= + (2k) 1) E a Tee | Retemecone (42) 
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The general case (no restrictions on v or V/v) is more difficult, but we can get 
an idea of the magnitude of R as follows. We can write 


ETAT 8 cos ddd 
to [$o2 + 2kV cos f + 2kv($ — bo)]" 
and for (18) to be valid the magnitude of R must be small compared with 26,2 as 


given by (18). By a similar procedure to that used above we can show that the 
average value of the integrand of (43) will be nearly zero, except over the first few 


phase revolutions (this is so even if ¢) = 0, and also for particles which pass through 
$= before being captured, since in this latter case, if V is increasing slowly, ¢ 
will be approximately an even function of ¢ near 6=0). Now, in this initial 


period, V will be small, and if bo is zero or positive (i.e. corresponding to a particle 
initially on or inside the equilibrium orbit), we can get a reasonable idea of the 
magnitude of R by writing 


R~2kV | 


RZORV | peace Cp aes 
‘ bo [bo + 2ko(p — d)]} 

since this integrand has nearly zero average value except over the initial part of 
the range. For ¢) negative we are still to some extent justified in accepting 
(44) since even though ¢ is decreasing we shall get approximately the same value 
over a few phase revolutions whether we integrate over a positive or a negative 
range of ¢, although in either case it must be realized that except when V,/v is 
fairly large we shall get, at best, only a rough estimate of R. 


If we accept (44), the magnitude of R becomes 
2kV 
l= Fe), hte ws 45 


where 


F(z) = {2x[(C( 00) — C(z))? + (S( 00) — S(z))?]}3 ; 


C and S are the Fresnel integrals and z = bo2/2kv. For large z, F(z) 27?, where 
the approximation is satisfactory for z>5. 
It is convenient to compute R in terms of a quantity 5 defined by 
kvt, =5(2RV,9,)*, 
Le. TA 2RV ea 200g ee yee: (46) 
which with (19) defines uy for the extreme particles captured. It is clear that 


8 is the ratio of the average orbit contraction during ¢, to the half-width of the 


region from which particles are captured. 
The condition that R be negligible is that |R| <4kV,®,, which, with (45), 
gives 
F(z) (Vi\3 ye 
tS a (— ] nate ees (47) 
t,(2RV,) > 20, (=) \ lp 
Further, for the extreme particles, 
B= (be 0) O1 142, 
the most unfavourable value being 
(ll 2) UN i a a eer (48) 
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Using (46) and (48), the value of t,(2kV,)? to make (47) an equality has been 
computed for a range cf values of V,/v and is plotted in Figure 8. The appro- 
priate condition that R be negligible is just that t,(2kV,)' has a large value 
compared with the ordinate of Figure 8 corresponding to the given value of 

TAGs 


10-0 


10 100 
Vi /v 
Figure 8. t(2kV,)? to make (47) equality. For the capture theory to be valid, t,(2kV,)* 
must be large compared with this value. 


These restrictions on ¢, leading to the capture theory developed in the text 
represent one extreme of the capture process, just as the instantaneous rise of 
the R.F. represents the other extreme; between the two there will be a continuous 
gradation between the capture of all particles within certain radius limits irrespec- 
tive of the initial phases, and the capture of only those particles which are initially 
within the synchrotron stable region. 
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Experiments on Electron Capture and Phase Stability in a 
14 Mev. Synchrotron 


Bye eRe SE ROANDS) oI, LUCK, 
Clarendon Laboratory, Oxford 


Communicated by Lord Cherwell ; M/S. received 11th Fuly 1949 


ABSTRACT. Experiments were made by interrupting the radio-frequency accelerating 
voltage for intervals during the acceleration in a synchrotron, and observing the effects on 
the output. The process of capture is found to be more efficient than would be expected 
for an instantaneously developed radio-frequency accelerating voltage. ‘The synchrotron 
acceleration may be interrupted for a short period without losing more particles than would 
be accounted for by contraction of the beam radius during the interruption, so that for 
sufficiently short periods of interruption the loss is negligible. "This demonstrates that the 
phase at which electrons enter the resonator during capture is of no consequence under the 
conditions prevailing. ‘The results are in quantitative agreement with a theory put forward 
by Kaiser which takes into account the finite rate of rise of resonator voltage, and with the 
general theory of synchrotron stability. 


UNGER OMe LOIN 
EASUREMENTS Of the efficiency of capture of electrons into the stable 
orbits of a synchrotron have been made using a 14 Mev. electron syn- 
chrotron (Abson and Holmes 1947). 

In this synchrotron the magnet is energized by 50c/s. a.c. with a D.c. bias 
current, the electrons being injected at an energy of about 10kv. and accelerated 
by betratron action to a total energy of about 2-5 Mev. when the betatron flux is 
beginning to saturate. At this point, with an equilibrium orbit radius of about 
7:-5cm., R.F. power is supplied to the resonator and synchrotron acceleration 
begins. The r.F. peak voltage for almost complete capture is found to be 130v. 
Other parameters which will be of interest are: 

Maximum voltage per turn required by the equilibrium 


particles . . - = i 3 fs as Up = 1-87 v. 
Magnetic field exponent .. oe kas cf ce n= 0°71 
Resonator frequency 5s aa EK a 2 Oem —O4+3- Neils: 
Equilibrium radius Pe “e = : t= 7) Cm. 
Width of orbit space at injection time .. sie 2m 5 cm. 
Resonator voltage Be : : e. V=30ive 


The notation of the previous paper (Kaiser 1950) is used throughout, and the 
paper referred to as I. 

The technique employed was to suppress the R.F. for a time 7 (variable from 
about 3 to 100 psec.) at various times T (corresponding to various values of the 
magnetic guide field H,) during the acceleration, and to examine the dependence 
of the x-ray output intensity at 2, = 14 Mev. (which is proportional to the number 
of particles in the beam) onr. The complete cycle of events was displayed on the 
screen of a cathode-ray tube, as in Figure 1, horizontal deflection being proportional 
to the magnet current. 

The x-ray intensity was measured by the amplitude of the pulse from a propor- 
tional counter placed in the beam, the counter having been calibrated against 
an ionization chamber type dosemeter. The dosemeter itself could not be used 

57-2 
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because of the x rays produced earlier in the cycle by the particles lost during the 
period of R.F. suppression. Since the master oscillator of the R.F. generator was 
maintained and only the power amplifier suppressed, there was phase coherence 
before and after suppression. Great pains were taken to ensure that there was no 
R.F. ‘break-through’ during the suppression period by suitably biasing the power 
amplifier, since it was found that quite a small break-through could give spurious 
results, in particular as regards the time, 7), required to suppress the beam 
completely. 

The shape and length of the negative blanking pulse were measured with a 
crystal rectifier loosely coupled to the resonator, and a cathode-ray oscilloscope 
containing a generator of 1 psec. timing pulses. The blanking pulse showing the 
timing markers is sketched in Figure 2. Both the decay and rise times were of the 


ee eel ae 


Figure 1. Sequence of operations during the acceleration cycle: a, betatron injection puise ; 
b, instant of establishing the R.F. ; c, X-ray output due to the few betatron electrons not 
captured ; d, R.F. suppression ; e, X-ray output due to electrons lost during suppression ; 
f, final R.F. suppression ; g, X-ray output due to 14 Mev. electrons. 


a ear 


Figure 2. Sketch of the negative blanking pulse with 1 usec. timing pulses superimposed. 


order of lsec. The actual instant of suppression was known from the calibration. 
of the variable delay in the blanking pulse generator (which was triggered at the 
instant of transition from betatron to synchrotron acceleration) and was related to 
the magnet current, and thus to H,, by measuring the displacement of the blanking 
pulse along the time-base of Figure 1. ‘This also enabled the r.F. duty cycle to be 
determined, and the rate of increase of H, to be checked for each experiment. 

The R.F. power input to the resonator was measured with a bolometer bridge 
and a wave selector coupled to the transmission line from the power amplifier. 
In all of the experiments described, this power was adjusted to give almost complete 
capture, in every case corresponding to a peak input of 28 watts. The measured O: 
of the quarter-wave resonator was 300 and its characteristic impedance 0-8 ohm, 
from which we can calculate for the shunt resistance a value of 306 ohms. Using 
the measured power input we find the peak resonator voltage V to be 130 volts. 
This value was expected to be very close to the actual peak voltage per turn 
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experienced by an electron, since the conducting surface of the resonator was 
continuous.* Immediately following the suppression we shall have a beam of 
electrons spread in radius and phase (relative to the r.F. phase) in a manner 

- dependent on the time of decay of the r.F. voltage. ‘The mean orbit radius will 
begin to contract at a rate 7, the beam keeping substantially the same width so long 
as the change in mean radius and the beam width are small compared withr,. The 
electron density (averaged with respect to 4) will remain substantially the same 
function of radius, but since particles on the inside and outside of the beam will 
lose and gain phase respectively, relative to those in the centre, the beam will rapidly 
become uniform in phase. On re-establishing the radiofrequency the region from 
which particles may be captured will overlap only a fraction of the beam, this fraction 
being captured and the remainder spiralling to the inside wall of the orbit vessel. 
We can examine the fraction recaptured as a function of 7 and compare it with that 
expected, both on the capture theory of paper I and on that for instantaneous R.F. 
decay andrise. Further, from a knowledge of + and 7), we can obtain a measure- 
ment of the width of the region from which particles may be captured at various 
times during the period of acceleration, enabling quantitative comparison with the 
value expected from theory. 


$2. INSTANTANEOUS FALL AND RISE OF THE R.F. VOLTAGE 
Immediately following suppression the bunch shape will be approximately 
given by 


eM 


=cos}¢, since VSv 


(where ¢, is the maximum half-width of the stable region prior to suppression) 
and is plotted in Figure 3(a). As the mean orbit radius decreases the bunch will 


Figure 3. 


(a) Bunch shape immediately after instantaneous suppression of R.F. voltage. 
(b) Bunch shape when the phase of the extreme electrons relative to the central ones has advanced or 
retarded by 5z radians. 


spread in phase in the manner indicated in Figure 3 (b) which shows the change in 
shape after the phase gain (or loss), 54, of the extreme electrons in the bunch 
relative to the central ones has reached 57 radians. 


* The resonator dimensions were such that there was no need to combat eddy currents by 
subdividing the conducting surfaces. 
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This phase gain can be shown to satisfy 


pier 


eM 


where V is the peak voltage prior to suppression and 67 is the amount of contraction 
of the mean orbit radius. We should thus expect the beam to become uniform in 
phase when 67/e,,;7v/2V, which is equivalent to 

ie i 

7) > AY 
since 67/e,;=27/7). For this machine, using the least favourable value of v/V, 
(1) becomes 7/7, > 0-01. 

When the radiofrequency is instantaneously re-established, the stable region 

will be of the usual form and will overlap a fraction of the beam. If we adopt the 
simplifying assumption that the beam becomes uniform in phase immediately after 
suppression (and that prior to this electrons were uniformly distributed throughout 
the stable region), the fraction of the electrons recaptured is easily calculated and is. 
plotted in curve (a) of Figure 4, giving the expected relative x-ray intensity Ifa s. 
as a function of 7/7). In fact, we would expect something approximating to the 
dotted curve of Figure 4, dropping rapidly from unity and approaching the curve (a) 
by a damped oscillation which becomes negligible when (1) is satisfied. 


1-0 


0 0-5 i) 
T/To 
Figure 4. Theoretical curves for relative x-ray output intensity against period of R.F, suppression. 


(a) Instantaneous suppression. 
(b) and (c) Finite r.F. decay and rise times. 


§3. FINITE RF. DECAY AND BUILD-UP TEMES 
The actual decay and build-up times ¢ were of the order of 1 wsec., and using 
the figures given above we find that t(2kV)* varies from about 30 to 15 over the 
period of synchrotron acceleration. In the same period V/v decreases from 


several hundred to about 70, so that equation (23) of I, 1.e. 2V/u>t(2kV)!> 4, is 
fairly well satisfied at all times. 
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In this case (on the theory of the previous paper) the bunch shape will be 
destroyed during the r.F. decay time and the half-width of the beam will contract 
to 1//2 times that of the stable region prior to suppression, i.e. to €y,/\/2. On 
re-establishing the radiofrequency all particles will be trapped which lie between 
the limits 7, + €,/+/2, so that we can evaluate the (intensity, 7) curve if we know the 
density—radius distribution of the electrons during the suppression (which may be 
similar to that in the original betatron beam). The curves (6) and (c) of Figure 4 
give the expected law for two possible electron distributions, the former on the 
assumption of a density function proportional to (1 —e?), where « is the displace- 
ment of the particle from the centre of the beam, and the latter on the (less likely) 
assumption that the density is constant throughout the beam, falling abruptly to 
zero at the edges. 

In the above, the effect of phase damping during the period between capture 
and R.F. suppression has been neglected, as has also the change in H, and w during 
the period of suppression, the latter being negligible during the short time interval 
involved. 


94. EXPERIMENTAL RESULTS 


A sample set of experimental results for ///,) plotted against 7/7, at various times 
during the acceleration is givenin Figure 5; this comprises three experimental runs 
taken on different days. The full curves are identical with (6) of Figure 4, while 
the broken curves represent an attempt to take into account the phase damping, 
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Figure 5. Experimental results for relative x-1ay output as a function of period of R.F. suppression, 
showing agreement with the theory. For convenience the scales of 7/7) for the various 
curves overlap, so that for each full curve 7/7) has the values 0 and 1 corresponding to 
ordinates 1 and 0 respectively. 


and are obtained from the full curves by assuming that the beam width during R.F. 
suppression is smaller than that determined from equation (24) of I by a factor 
proportional to H>*. 

The results show that the capture is nct of the form to be expected for an 
instantaneous R.F. decay and rise, and support the conclusion that the initial phase 
at which the particle enters the resonator is of no consequence, all particles within 
certain radius limits being captured. The agreement between the experimental 
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and theoretical results seems satisfactory, but is actually not as important as the 
general form of the results, since any phase dependence would lead, as 7 increases, 
to arapid fall of //J, from unity followed by a slower decline to zero (as, for instance, | 
curve (a) of Figure 4). 

The derivative of the experimental curve gives the form of the density—radius 
distribution in the beam during suppression, and an examination of the results for 
T = 160 psec. (and to a lesser extent 7 = 560 usec.) seems to indicate a slight density 
minimum near the centre of the beam. The fact that this does not appear in later 
curves may be due to space-charge effects and indicates that it probably derives 
from the original betatron beam. 

From a knowledge of ¢ during suppression, and of 7», we can calculate the radial 
width of the region from which particles may be captured and can thus obtain 
quantitative verification of the theory. 

From equations (7) and (24) of I, the expected half-width « of the beam during 
suppression is 


ib eV \3 
6 = oe ee (2) 
Bs2(1 —2) \wE,K 
where V, of course, is the voltage immediately prior to suppression. 
During this period the orbit radius of a particle in terms of its momentum, P, 1s 


pe 
T= He eee = Sw fe Brean: (3) 
where C is the velocity of light, so that 
Bate eee 
ee ek Ge 


Now f is not zero, due to betatron acceleration, and H is a function of 7, thus 


pe nee ee 
p. = Be E. and ie =f7. =nH, ry. 

where E is the rate of increase of energy due to betatron action and is related 
to v by 


2a Zak 
v= _—. 


Cc ew, 
From (3) he cals 
‘ $ Ws 
6p ; =} €Cv 
ivin ee ee 
giving r nb En), nee, ee (4) 
: ten: x: CCT VU 
and since «=7 97/2, we get «= ta(L n\n oe) a ae (5) 


Although we can evaluate « from (5) and compare it with the value obtained 


from (2), the following equivalent procedure has been found convenient. From 
(2) and (5) we get 


©. e| (e) «fe 


where 7 =4(7ceV/Kw,3)}, i.e. V= —— £:S.U; (7) 


Electron Capture and Phase Stability in a Synchrotron 73 
If » is measured in volt.usec.gauss~!, (7) becomes 
See CLM. ia” pe ee en are (8) 


In evaluating the results, the particles have been assumed extremely relativistic, 
i.e. K=1/(1 —n) throughout the acceleration. 

The full curve of Figure 6 gives v as a function of H, as determined by Abson 
and Holmes (1947)*, while the experimental points are those obtained by sub- 
stituting the measured values of ty in (6) and choosing 7 to give a good fit over the 
flat top of the (v, H,) curve (H,>2,500). The values of 7, were obtained by plotting 
I/TIy against 7, and for the first few curves the intersection with the horizontal axis 
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Figure 6. Experimental values of v, as a function of Hs, compared with the result obtained by 
direct measurement of the residual betatron induction and the rate of increase of Hs. 


was taken. ‘The later curves correspond to a period of the acceleration when v 
is almost constant, as also is ty; in order to compare the results with the curve of 
Figure 6 the relative values of 7) are more important than absolute values, so that, 
to obtain greater consistency it was assumed that the graphs of J/J) against 7/79 
were approximately anti-symmetrical about the line 7/7) =4, 7) thus being deter- 
mined from 


i ' (I) Iq) dr 


by graphical integration. The results so obtained agreed with the intersection of 
the curve with the horizontal axis to within the limits of estimation, but the 
consistency should be somewhat better. 

Using the values of 7 as determined above, the resonator voltage V was 
evaluated from (8) for each of the sets of results of Figure 6. The values are given 
in the Table, together with the value determined from the resonator power input. 


* A numerical error in Abson and Holmes’ paper of a factor 1/2 in the value of v has been 
corrected. 
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The agreement is well within experimental error, the scatter inthe magnitude of V 
determined from (8) corresponding to a scatter in 79 of less than | psec. 


Experiment (see Figure 6). . ° ~< A 

Nl <nets as oe 0-49 0-495 0-52 
V (volts) (from (8)).. a 122 124 136 
V (volts) (from power input) 130 130 130 


An extra check on the validity of the capture theory may be obtained by calcu- 
lating the radial zone of acceptance and comparing it with the expected width of 
the betraton beam at the instant of transition from betatron to synchrotron. For 
a resonator voltage of 130 v. the half-width of the acceptance zone is «y = 0-06 cm., 
while the assumption that the initial orbit space of width 5cm. is completely 
filled at injection, following which the beam contracts approximately as H~* 
(Kerst and Serber 1941), leads to «)=0-1cm. for the half-width of the betatron 
beam (neglecting free radial oscillations) immediately prior to transition. 


§5. CONCLUSIONS 

Experiments have been made on the release of electrons from, and their 
re-capture into, stable orbits in a synchrotron. The capture is found to be more 
efficient than would be expected by the application of the usual synchrotron 
equations for an instantaneously developed radio-frequency accelerating voltage. 
The results are in agreement with a theory (given in the previous paper) taking into 
account the finite rate at which the r.F. voltage is developed and with the general. 
theory of synchrotron stability on which the capture theory is based. In particular 
the form of the curves of Figure 5 indicates that the initial-phase at which electrons. 
enter the resonator is of no importance, while the agreement between the values of 
resonator voltage in the Table supports the conclusion that the width of the 
electron beam which can be captured is 1/,/2 times the width of the synchrotron 
stable region (it should be noted that the voltage V, calculated on the assumption 
of instantaneous capture, would be one half that given in the Table). The 
agreement between the experimental points and the curve of Figure 6 shows that the 
width of the stable region varies with magnetic field (or particle energy) in the 
manner required by the general theory. 

A modification of the experimental technique that we have used suggests itself 
for a direct measurement of the distribution of particles with radius in a betatron 
beam prior to transition to synchrotron acceleration, namely by measuring the 
output intensity as a function of the time of transition. Provided v remains 
constant over the interval concerned, the form of this distribution is just 
the derivative of the (x-ray intensity, time) curve between maximum and zero 


intensity, while if v does not remain constant, its variation, if known, can be 
taken into account. 
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The Electric Conductivity of Simple Semiconductors 


E. H. Putley, in this journal, has recently (1949) drawn attention to a series of papers in 
which Shifrin (1944) discussed the electric properties of impurity semiconductors, to be 
expected on the basis of the model originated by Wilson; Shifrin extended the limiting 
relations originally given so as to cover a much wider range of temperatures and impurity 
densities. In particular Shifrin has shown that for simple semiconductors the conductivity 
for the whole range of temperatures should rigorously be given by ; 


GO wlian (et A) ee, Wide en eae (1) 


where 4 is the quantity originally introduced by Sommerfeld (1928) and related to the 
thermodynamic potential € by In dA=y=-¢/RT. Further, he has shown that A and hence 
the conductivity has a maximum A), at a temperature Ty, given by 


Inge ARTS inAgkTo= 1 9 Fr: (2) 


‘These expressions refer to non-ionic lattices for which the mean free path of the electrons. 
is independent of their energy and inversely proportional to the temperature, with only one 
kind of impurity centre supplying the electron gas. They take fully into account its 
degeneracy. It has since been shown by a number of authors (Johnson and Lark-Horovitz 
1947, 1942-1949, Busch and Labhart 1946, Putley 1949) that this degeneracy plays in fact 
an important réle. 

The question then arises how to apply equations (1) and (2) in interpreting measured 
conductivities in the light of the relation 


n=np/{1+A exp (Ac/kT)}=2C(RT)*7fA(A), iw se ee (3) 
where 
2 x? dx 


gp Aa 
ny the number of impurity centres, 7 the number of free electrons, Ae the energy difference 


between the impurity level and the lower edge of the conduction band, 


C=4n(2m)3/2/h3 ~3-4 x 10° c.G.s. units. 


f(A) =3 


Tables of fy have been given by McDougall and Stoner (1939) but a tabulated form of fy 
does not lead to an explicit expression for A in terms of 7p, Ae and T. 
Now, Busch and Labhart have proposed to approximate f,; by an expression of the form 


fa=Art/41+bA+cA®), nee (5) 


valid for 4<5, i.e. in the range which one expects to find in semiconductors. In following 
this up, however, Busch and Labhart use instead of Shifrin’s formula (1) the expression 
we 
ax eae Sg 
which would only be justified if for semiconductors the time of relaxation 7 were independent 
of the energy of the electrons, but which is incompatible with the accepted conclusion that 
the mean free path is independent of the energy. 

Putley, on the other hand, uses Shifrin’s relation in connection with the tables of 
McDougall and Stoner to find 9, 7p, Ae. Inconsistently, however, he suggests obtaining oo. 
from the low-temperature range, whereas the low-temperature approximation of A leads 
directly only to Ae as the negative slope of In o(kT)? plotted against 1/2kT. From Ae and 
Tn, Am can be obtained from (2) and hence mp from (3) and oo from (1). But lacking an 
analytic expression for A, it is not easy to calculate o generally. 

Busch and Labhart, however, had obtained an analytic expression for A. Following 
up their suggestion one readily verifies that 


POO 2s Ay) ae eee ee (6) 
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reproduces the values of McDougall and Stoner for A< 6 with an accuracy better than 5%, 
which is quite satisfactory in the case of semiconductors. With (6), the temperature 
dependence of the conductivity of semiconductors becomes very clear. Equation (3) leads 
-on the one hand to 


nN=0°886C(R ID)? 2 /(O:25 2s) 5 nner were (7) 
and on the other hand to a quadratic equation for A or A™, 
A1=0-443C(RT)??/np—0-125 
+ [(0-443C(RT)3/2/ny,—0:125)?+-0-886C(RT)? exp (Ac/RT)/ny 3, 


so that, using (1), an explicit expression for o can be written down. 

When (8) is specialized for high or low T and the resulting value is inserted into (1) and 
(7), the well-known limiting expressions for the conductivity and the number of free electrons 
are obtained. Equation (2) follows at once from (8) by equating 0471/eT=0. Further 


Np =0-883C(kTm)?[1+Am exp (Ae/RTm)]/(0-25+Am). «----- (9) 


Hence, by combining an approximate expression for f; of the Busch and Labhart type 
with Shifrin’s results, an analytic expression for the conductivity of semiconductors valid 
for all temperatures becomes available, as well as a simple procedure for obtaining the 
constants. It should, however, be borne in mind that it is derived under the assumptions 
that the mean free path ot the electrons is determined exclusively by collisions with acoustic 
lattice vibrations, that only one kind of impurity centre supplies the electron gas and that 
conductivity takes place in one: band orly. 


Birkbeck College, W. EHRENBERG. 
University of London. 
1st November 1949. 
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REVIEWS OF BOOKS 


The Freezing of Supercooled Water, by N. Ernest Dorsey. Poe S2« 
Reprinted from Transactions of the American Philosophical Society, 38,. 
247, (Philadephia: American Philosophical Society, 1948.) $1.75. 


In 1938 Dr. Dorsey published a preliminary account of experiments on the supercooling 
of water. In this he showed that the extent to which a given sample could be supercooled 
without freezing was determined by the presence of foreign particles or nuclei, and that these 
could be removed or influenced by further purification, by allowing the sample to stand, or 
by heating. Very pure water could be supercooled to —21°. The determining factor was. 
thought to be the size of these nuclei. The present longer (82-page) account describes his 
further experiments carried out intermittently over the last ten years. The style is more 
discursive than is usual in an account of original research but he presents a large number of 
results which undoubtedly help to fill in the complicated picture which has arisen from his. 
earlier work and from the more recent work of Cwilong, Rau (extended to much lower 
temperatures), Weickmann, etc. Dorsey’s experiments are all measurements of top, the 
temperature of spontaneous freezing, and for water in his temperature range this means that 
they are investigations of nucleation, the rate of crystallization being always very large after’ 
initiation. He has investigated carefully the effect of prolonged maintenance at particular: 
(higher) temperatures and of various other forms of pre-treatment, and has also established the 
existence of ‘ preferred temperatures ’ for the freezing of his samples. Only speculative 
suggestions are offered about the actual nature of the nuclei, and some of his experiments are 
of qualitative rather than quantitative interest—as for instance his tests on “‘ ordinary 
distilled water containing an unknown amount of alcohol and an unpervforated glass bead ’” 
and “‘ vacuum distilled water into which ethy! alcohol from the bath leaked through a crack’’.. 

While the experimental results presented are useful and extensive, the theoretical 
discussion is of less value. Dr. Dorsey is critical of previous theories of nucleation, par- 
ticularly of the ‘homogeneous’ theory arising from the work of Gibks, J. J. Thomson, 
Frenkel, Becker, etc., which considers the onset of crystallization to be due to spontaneous. 
nucleation. Favouring in general some form of ‘ heterogeneous’ theory, that is, one in. 
which nuclei already present are assumed to be responsible for crystallization, he has (to use 
his own phrase) “ worked up’ a new detailed theory in which ice embryos are assumed to be 
formed in various ways from water molecules adsorbed on larger foreign paiticles. He is 
probably right in saying that the ‘ homogeneous ’ theory is inapplicable to his experimental 
results, but is unduly polemical on this and other points. G. O. JONES. 


Foundations of Modern Physics, by T. B. Brown. Pp. xvi+391. Second 
Edition. (New York: John Wiley and Sons; British Agents: Chapman 
and Hall, 1949.) 30s. 


Professor Brown’s book is successful if its aim is to give to science students, in their last 
year at school or in their first year at the university, a sketch of the type of experiments and 
an outline of the ideas which form the basis of modern physics. 

The production of the present work is excellent, the printing is good, the diagrams 
are clear and concentrate on fundamentals, and there are extremely few errors of any 
kind. There does not seem to be any way in which the treatment of so wide a subject 
could be substantially improved within the limits of 400 pages. Yet it does not seem 
possible to deny that the analysis is shallow ; also the contribution of American workers 
is over-emphasized in the last few chapters. Perhaps the author can hardly be blamed 
for this since the advance of modern physics has been so rapid and on so wide a front 
that some selection is inevitable. Indeed, the reviewer feels that the days when a semi- 
popular exposition of physics can be undertaken are rapidly drawing toa close. No one 
expects advances in modern mathematics to be made generally intelligible or available 
to first-year students, and certainly the number and quality of those engaged in advancing 
modern physics is such that the breadth, content and ccmplexity of the subject are likely 
to increase still further. F. C. CHAMPION. 
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Einfiihrung in die Atomphysik, by W. FINKELNBURG. Pp. vii +339. (Berlin : 
Springer-Verlag, 1948.) 28 D.M. 


The appearance of this excellent text by Professor Finkelnburg is a welcome sign of the 
revival of scientific activity in post-war Germany. The book is based on courses given for 
ten years to classes of physicists, chemists, engineers and biologists. It is approximately 
of the standard required for a first degree in physics in this country. ‘There are now a 
number of works in English covering similar subject-imatter of about the same standard, but 
this book compares favourably with any of them in breadth and interest of presentation. — 

Professor Finkelnburg aims at presenting atomic physics, not as a series of isolated topics, 
but as a subject with an underlying internal unity. Chapters of the book deal with atomic 
structure and atomic spectra, nuclear physics, molecular physics and the physics of the solid 
and liquid states. The outlook is predominantly experimental, and a chapter on the basic 
ideas of quantum theory enable the author to show how the multitude of experimental data 
-on atomic physics, which twenty-five years ago seemed to defy attempts at rationalization, 
can now be reduced to a rational form in terms of quantum mechanics. : 

The book as a textbook is thoroughly up to date and contains sections on nuclear fission, 
‘isotope separation and the betratron. 

It is attractively printed and shows that the great reputation built up by Messrs. Springer 
as publishers of scientific literature in former days is likely to be maintained in these post-war 
‘years. E. H. S. BURHOP. 
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The Elastic Constants of a Solid containing Spherical Holes, by J. K. MackEnzir. 


ABSTRACT. The effective bulk and shear moduli are calculated by a self-consistent 
method due to Fréhlich and Sack. The bulk modulus k is determined by applying a 
hydrostatic pressure, and the shear modulus jz by applying a simple homogeneous shear 
stress, to a large sphere. Each hole is surrounded by a spherical shell of real material, and 
the reaction of the rest of the material is estimated by replacing it by equivalent homogeneous 
material. For consistency, both the density and the displacement of the outer spherical 
boundary must be the same whether the hole and its surrounding shell are replaced by 
equivalent material or not. The effective elastic constants calculated from these conditions 
are 
1/k=1/kop +3(1—p)/4u40p + O[(1 — p)*1, 
(Ho —P)/Mo = 5(1 — p)(3ko + 4p0)/(9Ro + 80) + OL — p)?], 


where kp and [uy refer to the real material and p is the density of the actual material relative 
to that of the real material ; in the next approximation k depends on the standard deviation 
of the volumes of the holes. 

The dilatation due to a distribution of pressures in the holes is p(1/k—1/ky), where p 
is the mean obtained when the pressure in each hole has a weight proportional to the volume 
of the hole. By using the hydrodynamic analogue of the elastic problem, the theory is 
briefly applied to the theory of sintering, and used to discuss the effective viscosity of a 
liquid containing smali air bubbles. 


On Movements of Small Ferromagnetic Particles in Inhomogeneous Magnetic Fields, 
by FRIEDRICH BLAHA. 


ABSTRACT. Ferromagnetic particles of diameter from 107° to 10-> cm. when exposed 
to light move along the lines of force of an inhomogeneous magnetic field. It does not seem 
possible to explain the movement when only dipole properties of the particles are taken into 
account. Field strengths of 10~? gauss are sufficient to cause distinct movement. ‘The 
motions resemble those of electrically charged particles in electric fields. 


Dead Times of Self-Quenching Counters, by B. COLLINGE. 


ABSTRACT. Experiments on Geiger counters are described in which localization of the 
ion sheath was brought about by reducing the wire potential immediately after each count. 
The effect on the dead time is discussed. The dead time of a counter 30 cm. long was 
found to be 20 psec. and independent of the counting rate up to 1°8 x 10® counts per minute. 
‘The reduction of the dead time of short counters was also investigated. 


Breakdown in Cold-Cathode Tubes at Low Pressure, by F. G. HEYMANN. 


ABSTRACT. When, in cold-cathode discharge tubes, electrode spacings and surface 
areas deviate from the ideal of infinite parallel planes, it is found that the breakdown voltage 
between electrodes increases. This is explained as a decrease in the effective 7 of the gas 
due to loss of electrons and positive ions by diffusion to the walls of the container. ‘The loss 
factor per unit potential difference is shown theoretically to be inversely proportional to the 
field strength and tube radius, although this is not fully verified by experiment. 

Paschen curves obtained experimentally for potassium and nickel cathodes in argon and 
in a neon—argon mixture at low pressures are shown and from these, values of y for the two 
cathode surfaces are obtained as a function of E/p). The apparatus used for measuring 
breakdown voltage is described. 

Earlier theories of statistical and formative time delays are extended to cover the case of a 
rising overvoltage and also the case where the primary electrons appear in bursts as with 
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ionization by a-particles. The shape of the statistical distribution curve is an indicaticn of 
whether the primary electrons have been produced singly or in bursts. The overvoltage 
AV of breakdown due to formative lag and the rate of rise of the uniformly increasing 
applied voltage V are found to bear the relation AV “1/(dV/dt). This has been experi-: 
mentally verified. 


The Behaviour of Multiple Circuit Magnetrons in the Neighbourhood of the Critical 
Anode Voltage, by W. E. WILLsHaw and R. G. ROBERTSHAW. 


ABSTRACT. It is pointed out that the mechanism of operation of the multiple circuit 
magnetron oscillator in the region of minimum magnetic field and voltage, where the 
efficiency is commonly assumed to approach zero, should approach that of an oscillator of 
the travelling wave tube type, providing that a cathode of suitable size is used. Useful 
efficiencies should thus be obtainable under these conditions. 

Details of experiments are given in which an electronic efficiency of 12% was obtained at a 
wavelength of 3 cm. at values of magnetic field and voltage several times lowe: than those used _- 
for high efficiency operation. The mode of operation was determined by the value of the 
magnetic field, a given mode being maintained over a range of magnetic field of the order of 
8°4. The anode voltage was about 70°% of the critical value. 

The experimental results generally support the hypothesis and suggest that the minimum 
voltage regime should be of extreme importance for work at the highest radio frequencies. 


Measurements of the Reflection Coefficient of Water at a Wavelength of 8-7 mm... 
by-D.7G-KIECY, 


ABSTRACT. An account is given of measurements (carried out in March 1947) of the 
reflection coefficient of water at a wavelength of 8-7 mm. over a range of angles of incidence.. 

The method employed is to measure the relative field strengths of the direct waves and 
waves reflected from a trough of water using free-space propagation and high-gain aerials.. 

The following electrical constants of water have been computed from the measured 
results (water temperature 11:1°c.): refractive index 4:-40+0-24, dielectric constant 
10°86+2-21, absorption coefficient 2-91+0-06, conductivity/frequency 12-82+0-42, 
Brewster angle 79°. 


Coupling of the Ordinary and Extraordinary Rays in the Ionosphere, by 
T. L. ECKERSLEY. 


ABSTRACT. In this paper an approximate solution to the wave equation is given for 
propagation in an ionosphere in which the gradient of the density N is in the vertical, z, 
direction only, and in which account is taken of the earth’s magnetic field. It corresponds. 
exactly to the ray theory and expresses a quantity Z, which is the z derivative of the phase 
function S, by a quartic equation. Z can be represented as a function of ¢ (which is pro- 
portional to N) on a four-sheeted Riemann surface, and the branch points are studied for 
the case of vertical incidence for which Z becomes the refractive index. By considering 
the branch points in the complex ¢ plane, the amount of the coupling between the ordinary 
and the isolated extraordinary branches of the (Z, ¢) curves can be expressed as a function 
of the obliquity of the magnetic field. The triple splitting of rays reflected from the 
ionosphere, observable where the field is nearly vertical, can thus be explained, and the 
theory is substantiated by the observation that the polarizations of the echoes on the (P’, f) 
records are ordinary, ordinary and extraordinary in order of increasing critical frequencies, 
as given by the branches of the (Z, C) curves. 
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